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II 


A  small-perturbation  stability  analysis  of  a  doubly  infinite 
array  of  interdigitated ,  right  circular  helical  vortices  has 
been  formulated.  This  array  corresponds  to  the  vortices 
trailed  from  the  tips  of  the  blades  of  a  helicopter  rotor  or 
propeller  in  static  thrust  or  axial  flight  condition  and  at 
great  distance  from  the  plane  of  rotation  of  the  blades.  The 
analysis  makes  use  of  the  Biot-Savart  law  of  induction  and  the 
Vorticity  Transport  Theorem.  The  singularities  in  the  Biot- 
Savart  integration  for  self-induction  have  been  eliminated  by 
substituting  appropriate  approximate  functions.  Near- singular 
behavior  in  ether  integrals  has  been  minimized  by  adding  and 
subtracting  functions  with  similar  near- singular  behavior  and 
which  have  exact,  closed-form  integrals.  The  calculations  of 
induced  perturbation  velocities  and  those  required  for  the 
stability  analysis  have  been  programmed  for  digital  computer. 


Numerical  results  have  been  obtained  for  two-,  three-,  four-, 
five-,  six-helix  arrays  representing  the  vortices  trailing 
from  the  same  number  of  blades.  The  special  case  of  a  single 
helix  has  been  run  and  the  results  compared  with  those  present¬ 
ly  available  from  studies  by  Levy  and  Forsdyke  and  by  Widnall. 

A  continuum  of  instability  modes  has  been  found  associated 
with  all  values  of  wave  numbers;  only  modes  with  wave  number  0 
and  1  are  so  much  as  neutrally  stable,  and  only  for  the  case  of 
a  single  helix.  The  most  unstable  modes  involve  the  most 
axial  motion  of  adjacent  vortex  segments  relative  to  each 
other.  By  "adjacent  segments"  is  meant  vortex  segments  above 
and  below  each  other  (i.e.,  at  the  same  azimuthal  location)  on 
adjacent  coils  of  the  same  or  neighboring  helices.  Further¬ 
more,  the  larger  the  percentage  of  the  helical  arc  length 
involved  in  such  motion,  the  more  rapidly  the  distortion  will 
diverge.  Maximum  divergence  rates  in  the  unstable  modes 
increase  as  the  helix  pitch  decreases,  increase  as  the  number 
of  helices  increase, and  decrease  as  the  number  of  cycles  of 
deformations  in  one  turn  of  the  helix  (i.e.,  wave  number) 
increases.  The  larger  the  helix  filament  core  diameter,  the 
more  sensitive  the  analysis  is  to  the  means  by  which  the  singu¬ 
larities  in  the  self-induction  integrals  are  eliminated. 
Increasing  core  diameters,  however,  reduces  the  maximum  diver¬ 
gence  rates. 
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CHAPTER  1 


INTRODUCTION 


1.1  IMPORTANCE  OF  AN  ACCURATE  KNOWLEDGE  OF  WAKE  GEOMETRY 


The  geometry  of  the  wake  beneath  a  helicopter  rotor  is  impor¬ 
tant  to  the  accurate  prediction  of  the  airloads  acting  on  the 
blades.  Determination  of  such  airloads  is  prerequisite  to 
nearly  every  area  of  rotary-wing  design  including  performance, 
vibrations,  acoustics  and  structural  integrity.  An  accurate 
knowledge  of  the  wake-induced  velocity  field  is  also  impor¬ 
tant  for  determining  (1)  interference  effects  between  rotors, 

(2)  positioning  and  effectiveness  of  auxiliary  surfaces, 

(3)  interference  with  the  fuselage  and  other  non-lifting 
bodies,  etc. 

The  impact  of  wake  geometry  is  felt  directly  in  attempts  to 
predict  wake-induced  velocity  fields.  Inability  to  do  so  with 
acceptable  accuracy  has  limited  the  state  of  the  art  of  heli¬ 
copter  rotor  designs  for  perhaps  a  decade.  This  is  a  problem 
on  which  considerable  emphasis  has  been  and  is  still  being 
placed.  Some  of  the  reasons  for  this  emphasis  include  the 
following:  (1)  The  desirability  of  attaining  peak  lift  system 
performance  is  greater  with  helicopters  than  with  conventional 
aircraft  because  of  generally  less  favorable  payload  and  power 
to  gross  weight  ratios.  (2)  The  distinctive  acoustic  signature 
of  a  lifting  rotor  can  be  troublesome,  and  this  imposes  re¬ 
quirements  on  the  induced  velocity  field  that  have  no  counter¬ 
parts  in  fixed-wing  aircraft  design.  (3)  Not  only  is  the 
geometry  of  the  rotor's  vortex  system  considerably  more  compli¬ 
cated  than  that  of  a  fixed-wing  aircraft,  but  it  is  of  greater 
influence,  since  rotation  causes  the  wake  vorticity  to  remain 
longer  in  the  vicinity  of  the  blades.  (4)  Unsteady  effects 
must  almost  always  be  considered  since  velocity  at  the  blade 
element  varies  with  time  even  when  the  helicopter  is  in  steady 
flight.  (5)  Aeroelastic  effects  are  always  a  factor  since 
rotor  blades  are  relatively  more  flexible  than  conventional 
fixed  wings. 


1.2  CLASSICAL  METHODS  OF  COMPUTING  WAKE -INDUCED  FIELDS 


The  simplest  hovering  theory  is  the  so-called  actuator  disc 
analysis,  where  the  induced  velocities  are  calculated  in  terms 
of  momentum  and  energy  changes  as  the  flow  crosses  an  "actuator 
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disc".  See,  for  example,  Glauert  [Reference  1],  Shapiro  [Ref¬ 
erence  2]  and  Gessow  and  Meyers  [Reference  3],  Refinements 
of  this  approach  accounting  for  nonuniformities  around  the 
azimuth  have  been  made  by  Mangier  [Reference  4]  and  Loewy  and 
Joglekar  [Reference  5] .  The  next  step  taken  in  developing  an 
analytical  model  of  the  hovering  rotor  was  to  consider  the 
individual  blade  elements.  The  rotor  disc  was  divided  into 
elemental  annuli  through  each  of  which  the  momentum  change  was 
equated  to  the  blade  element  lift  for  a  given  blade  pitch. 

This  method  was  first  introduced  by  S.  Goldstein  [Reference  6], 
It  has  most  recently  been  developed  ♦o  account  for  experimental 
airfoil  characteristics  and  labeled  nonlinear  strip  theory" 
by  Jenny,  Olson  and  Landgrebe  [Reference  7], 


1.3  VORTEX  TUBE  REPRESENTATION  OF  HELICAL  WAKE  AND 

distorted  Tv'a'kE'  geometries  ~ 


Considerable  research  has  been  and  is  being  done  to  develop  a 
rotary-wing  equivalent  to  the  classical  lifting-line  type  of 
inflow  analysis  employed  with  success  for  fixed-wing  aircraft. 
Among  advances  made  in  the  last  decade  are  methods  of  wake 
analysis  which  deal  with  the  three-dimensional  array  of  vortic- 
ity  in  the  wake  "shed  and  trailed"*  from  the  finite  number  of 
blades  in  a  lifting  rotor.  Such  analyses  [Miller,  Reference  8; 
DuWaldt,  Reference  9,  etc.]  have  the  following  classical  bases: 
(1)  Making  use  of  Prandtl's  theory  of  airfoils,  the  wing  is 
represented  by  a  "lifting  line"  with  a  circulating  flow  around 
it.  This  circulating  flow  in  turn  gives  rise  to  circulation, 
r  ,  which  is  related  to  both  the  aerodynamic  forces  and  the 
strength  of  a  vortex  called  the  "bound  vortex"  which  is, 
essentially,  the  lifting  line.  (2)  The  aerodynamic  forces  are 
also  related  to  the  induced  velocities  at  the  lifting  line 
through  the  angle  of  attack,  on  which  dimensional  analysis 
shows  lift  to  be  dependent.  (3)  A  fundamental  law  of  the 
mechanics  of  perfect  fluids  (Kelvin's  theorem  or  the  Helmholtz 
theorem)  is  involved;  namely,  that  a  vortex  can  not  end  in  the 
middle  of  a  fluid,  so  the  bound  vortex  continues  in  the  fluid 
beyond  the  extent  of  the  blade,  as  a  "free"  vortex  trailing 
from  the  wing  tips.  (4)  Experience  has  shown  that  when  a 
sheet  of  free  vortices  trails  in  a  distributed  manner  from  some 
significant  spanwise  length  of  a  blade's  trailing  edge,  it 
then  rolls  up  into  a  pair  of  vortex  tubes  of  equal  strength  and 
opposite  directions  and  is  "transported"  away  from  the  blade 
which  trailed  it  by  the  rotational  velocity  and  the  downward 


*  "Shed"  vorticity  arises  from  time-wise  variations  in  lift, 
"trailed"  vorticity  from  spanwise  variations.  Only  the 
latter  kind  is  considered  in  this  report. 
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axial  induced  velocity  generated  by  t lie  rotor.  (5)  The 
velocity  induced  at  the  airfoil  by  both  bound  and  trailing 
vorticity  is  given  by  the  Biot-Savart  law  of  induction. 

The  classical  wake  is  generally  considered  in  two  distinct 
parts:  near  wake  and  far  wake.  The  near  wake  is  defined  as  a 
set  of  circular  helices  of  finite  length  trailing  from  the 
blade  tips  and  near  the  blade  roots,  together  with  the 
"lifting  line"  or  bound  vortex  representation  of  the  rotor 
blades.  The  length  of  this  set  of  trailing  vortices  is 
usually  in  terms  of  some  number  of  rotations  below  the  "lifting 
lines".  The  far  wake,  on  the  other  hand,  can  be  treated  as 
consisting  of  an  array  of  infinitely  long,  interdigitated  cir¬ 
cular  helices  as  generated  by  the  blade  tip  only. 

The  wake-induced  field  of  a  hovering  rotor,  therefore,  can  be 
inferred  from  proper  consideration  of  trailing  vorticity,  i.e., 
a  set  of  interdigitated  circular  vortex  helices  generated  by 
the  tip  of  each  blade  in  the  rotor.  The  question  is,  "How  can 
this  vortex  helix  representation  be  used  to  deal  more  funda¬ 
mentally  with  rotor -induced  flows?" 

Although  most  of  the  recent  methods  of  blade  airload  analysis 
represent  a  considerable  improvement  in  the  state  of  the  art, 
they  still  do  not  permit  accurate  calculation  of  blade  airloads 
and  associated  bending  moments.  The  fact  that  wake  geometry  is 
prescribed  in  advance  imposes  a  fundamental  limit  to  such 
methods.  Optimum  rotor  design  is  not  likely  to  be  achieved 
until  a  distorted  wake  geometry  can  be  accounted  for  rather 
than  prescribing  a  rigid  helical  geometry.  Tarrarine  [Refer¬ 
ence  10],  Landgrebe  [Reference  11]  and  others  have  shown, 
using  flow  visualization  techniques,  that  the  actual  rotor  wake 
geometry  differs  significantly  from  the  pure  helix.  Relatively 
recent  efforts,  therefore,  have  concentrated  on  the  development 
of  analytical  methods  for  predicting  the  distorted  helical 
geometry  of  the  actual  rotor  wake.  In  one  approach  [Reference 
5],  the  nonuniform  velocity  fields  predicted  by  momentum  theory 
have  been  used  to  calculate  the  distortion  of  a  canted,  but 
otherwise  pure,  helix.  In  others  [References  12,  13],  the 
self- induced  distortions  of  a  rotor  wake  represented  by  numer¬ 
ous  discrete  vortex  elements  are  computed  by  application  of  the 
classical  Biot-Savart  law  involving  both  numerical  integration 
and  iteration  techniques. 


1.4  MOTIVATION  FOR  THF.  PRESENT  WORK  AND  DESCRIPTION  OF  THE 

Senegal  approach  adoPtED 


The  determination  and  use  of  a  distorted  wake  geometry,  free 


from  the  constraint  of  conforming  with  a  pure  helix,  would  be  \ 
expected  to  provide  a  more  accurate  means  for  computing  the 
instantaneous  rotor  flow  field  and  the  associated  blade  air¬ 
loads.  It  seems  significant  to  this  objective  that,  in  tests 
and  wake  geometry  analyses,  evidence  of  apparent  instabilities 
in  the  tip  vortex  patterns  was  discerned;  in  experiments,  rapid 
deviations  from  the  helical  geometry  were  observed  [Reference 
10],  and  during  numerical  integrations, the  numbers  seemed  to 
diverge  in  successive  iterations  [Reference  11].  In  either 
case,  distorted  rotor  wakes  were  obtained.  In  Reference  11  it 
was  noted  that  the  portion  of  the  wake  in  the  immediate  vicini¬ 
ty  of  the  rotor  plane  (extremely  near  wake)  did  converge  in  the 
calculations,  with  the  degree  of  convergence  improving  with  in¬ 
creasing  proximity  to  the  rotor.  This,  presumably,  is  due  to 
the  fact  that  one  end  of  the  wake  is  tied  to  a  blade  whose 
position  is  prescribed.  Close  examination  of  the  computed  re¬ 
sults,  however,  indicated  that  the  wake  did  become  unstable  at 
moderate  distances  from  the  rotor.  It  therefore  appears  that 
to  gain  insight  as  to  the  completeness  needed  for  acceptable 
theoretical  rotor  blade  analysis  methods,  stability  of  the  wake 
should  be  studied.  Further,  it  seems  prudent  to  examine  the 
far  wake  first,  both  because  of  its  greater  simplicity  and  be¬ 
cause  of  its  apparently  greater  tendencies  toward  instability. 

Two-dimensional  flow  fields  behind  a  circular  cylinder  were 
found  by  Kfrmtfn  [Reference  14]  to  be  stable  only  when  vortices 
were  in  the  staggered  formation,  which  has  come  to  be  known  as 
the  "Kfrmtfn  Vortex  Street",  A  "Kirmin  Vortex  Street"  might  be 
thought  of  as  the  cross  section  of  a  right  circular  helical 
vortex  on  a  plane  through  its  axis.  Such  an  analogy  suggested 
that  the  interdigitated  vortex  helix  system  associated  with  two 
or  more  blades  would,  by  extension  of  Karman's  analysis,  be 
more  unstable  than  a  single  helical  vortex.  A  literature  sur¬ 
vey  revealed  no  analyses  of  the  stability  of  an  interdigitated 
vortex  helix  array  such  as  might  represent  the  wake  of  a  multi- 
bladed  rotor.  An  early  stability  analysis  of  a  single  helical 
vortex,  however,  was  carried  out  by  Levy  and  Forsdyke  [Refer¬ 
ence  IS],  This  work  was  prompted  by  flow  visualization  experi¬ 
ments  in  the  wind  tunnel  with  a  stationary,  circular  disc  with 
its  polar  axis  parallel  to  the  flow.  At  high  values  of  Rey¬ 
nolds  numbers,  it  was  found  that  its  wake  consisted  of  a  ro¬ 
tating  helical  vortex  filament. 

The  Levy  and  Forsdyke  analysis,  while  basically  like  that  of 
Karman,  involved  numerical  techniques  which,  without  modern 
computing,  led  to  some  erroneous  results.  Nevertheless,  the 
approach  adopted  here  for  the  stability  analysis  of  interdigi¬ 
tated  helical  vortex  systems  follows  the  same  general  line.  A 
system  of  helical  vortices  is  postulated.  Perturbation  dis¬ 
placements  are  assumed,  and  the  self-induced  and  mutual- induced 
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perturbation  velocities  at  the  perturbed  vortices  are  calcu¬ 
lated.  Subsequently  their  growth  rate  is  determined  using  the 
vorticity  transport  theorem.  As  was  found  in  an  independent 
study  of  the  stability  of  a  single  helical  vortex,  published 
after  the  initiation  of  the  research  reported  here  TNidnall, 
Reference  16],  the  vortex  system  is  unstable  under  most  situa¬ 
tions.  The  results  of  the  analysis  reported  here,  however, 
show  that  stability  and  growth  rate  of  perturbations  depend 
very  strongly  on  the  number  of  blades  in  the  multi-bladed  case. 


5 


CHAPTER  2 


THEORETICAL  ANALYSIS 


2.1  BACKGROUND 


Where  instabilities  exist,  iterative  techniques  will  not 
converge,  and  direct  numerical  integrations  with  respect  to 
time  will  be  greatly  enhanced  by  a  prior  knowledge  of  the 
unstable  inodes  and  their  rates  of  divergence.  As  noted  in  the 
preceding  chapter,  Levy  and  Forsdyke  [Reference  15]  considered 
a  single,  doubly  infinite,  constant-diameter,  helical  vortex. 

They  found  it  to  be  unstable  for  disturbance  modes  with  no 
distortion  for  vortex  helix  pitch  angles  less  than  0.3  radian. 
Their  analysis,  however,  suffered  from  several  shortcomings. 

The  numerical  integration  was  performed  using  a  planimeter, 
which  is  understandable,  considering  the  early  date  of  their 
work  (1928),  but  nevertheless  the  loss  in  accuracy  contributed 
to  errors  in  their  results,  and  the  approximations  used  in 
evaluating  the  singular  integrals  were  quite  crude. 

Widnall [Reference  16]  recently  investigated  the  stability  of  a 
single  helical  vortex  filament  with  a  finite  core, using  the 
method  of  matched  asymptotic  expansions.  The  results  of  that 
work,  except  for  certain  differences  in  evaluating  the  singular 
integrals  encountered  in  the  analysis,  were  obtained  as  a 
special  case  in  the  present  investigation.  In  another  relative¬ 
ly  recent  study,  S.  C.  Crow  [Reference  17]  dealt  with  the 
stability  of  a  pair  of  linear,  parallel  infinite  vortices  of 
constant  strength,  as  trailed  from  the  tips  of  the  wings  of 
fixed-wing  aircraft.  In  theory,  such  a  pair  of  trailing 
vortices  could  be  given  as  a  special  case  of  the  family  of 
interdigitated,  helical  vortices  considered  here;  specifically, 
a  set  of  two  interdigitated  helical  vortices  as  trailed  from 
the  tips  of  a  two-bladed  rotor,  but  with  infinite  pitch.  In 
practice,  the  approximations  used  in  performing  the  numerical 
integrations  over  space  precluded  such  a  quantitative  comparison. 
However,  Crow's  study  did  provide  insight  into  ways  of  looking 
into  the  mutual  inductance  effects  of  the  two  vortices  in  the 
present  investigation. 


2.2  CAPSULE  STATEMENT  OF  THE  WORK  DONE 

The  present  work  deals  with  the  stability  of  interdigitated 
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helical  vortices,  representing  the  part  of  the  far  wake  trailed 
from  the  tips  of  multibladed  helicopter  rotors  in  hovering  or 
vertical  flight.  A  review  of  the  literature  revealed  that  only 
References  IS  and  16,  which  are  limited  to  single  helices 
contained  pertinent  work.  The  analysis  in  this  report  involves 
the  following  steps:  (1)  An  n-bladed,  hovering  rotor  is  repre¬ 
sented  by  interdigitated,  doubly  infinite,  right  circular 
helical  vortices  of  constant  diameter.  (2)  The  system  of  n 
vortices  is  perturbed  from  its  normal  helical  configuration  and 
position  in  radial,  circumferential, and  axial  directions. 

(3)  Application  of  the  Biot-Savart  law  and  the  requirement 
that  the  time  rate  of  perturbation  displacements  be  compatible 
with  the  resulting  induced  velocities  yields  a  system  of 
eigenvalue  equations  for  the  admissible  characteristics  of 
these  perturbations. 


2.3  DEVELOPMENT 


2.3.1 Mathematical  Definition  of  Infinite  Helical  Vortices 


Referring  to  Figure  1 ,  which  makes  use  of  a  left-handed 
Cartesian  coordinate  system,  the  parametric  equation  of  a  cir¬ 
cular  helix  of  radius  r  and  pitch  k  is 


x*vcose  ,  ^  *  r  sin  e  ,  2:“  fere  >  -ooce<.oo 

To  mathematically  define  n  coaxial  vortices  of  radius  r  and 
pitch  k,  one  can  define  an  azimuthal  separation  angle  between 
the  first  helix  and  successive  helices  (see  Figure  1)  as 
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If  e_  represents  the  parameter  for  the  mth  vortex,  the  para¬ 
metric  equation  for  the  m*h  vortex  can  then  be  written  as 


Xyv>  -  r  COS  C9rv>  *■ 

m  r  s»n  (e™  4-  4U) 
*  ,v>  -  k  r  ©m 


2.3.2  Ctlculations  of  the  Induced  Velocities 


The  kinematic  relation  between  vorticity  and  velocity  in  an 
incompressible  fluid  is  given  by  the  well-known  Biot-Savart 
law.  This  equation  can  be  applied  to  the  system  of  n  vortices 
and  results  in  the  expression 


Tit  -  y  £ja  I  R*v,)>xdui  m 

4tt  \* — 

where  Itp  ■  induced  velocity  on  p1*  vortex. 

rm-  circulation  strength  of  rnh  vortex. 

elemental  length  vector  on  m’1'  vortex. 

F?«  relative  position  vector  of  point  on  p1K  vortex 
relative  to  points  on  m?' vortex. 

Referring  to  Figure  1,  a  vector  representing  an  element  of 
length  along  the  rrt*  vortex  can  be  written  as 


and 

Rwi^  m  (W Kp )  >  c^M-yp)  t 

Note  that  the  integration  in  equation  (1)  must  be  performed 
for  all  values  of  m ,  including  the  case  where  m»p,  to 
account  for  the  fact  that  the  p’’*'  vortex  induces  velocity  on 
itself.  Primes  have,  therefore,  been  used  to  identify  the 
general  points  involved  in  the  space  integration,  as  contrasted 
to  the  specific  point  at  which  the  induced  velocities  are 
evaluated. 

Rewriting  the  Cartesian  coordinate  description  for  the  m-"’  and 
pih  vortex, 

*-m  -  rcoa 
y*  -  r  sin 

«  ftre* 
and 

xp  »  r  Co*  t(?p  +  q>  ) 
yp=  r  Sin  COpf  y>) 

4P=  Rr©p 
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Therefore,  dL^  ®  [- r  Sm  (6m  +4m)  »  ^wi+^wO  ,  krj  de'm 

Rmp  -  L'C,rc<5St®^+MU)-)rc<5«sC0p+^)}  ,  {rsin(.Gm+^-  vsint,©^} 

t  \fcrte^-eP^J 

In  deterrainantal  form,  ^  x  d£m  can  be  written  as 


,  i-rcosC®^,  +  4>«)  {rsiACei«-t|w*)  *r<a^-qp) 

f  -rcos  CSp-Mfp)}  -  rsiyvC©p-t-»|V>} 


-  r  s.nte,;+4u) 


vcosce^+4U)  hr 


where  e  ,  e  ,  e  are  unit  vectors  in  the  Cartesian  x,  y  and  z 
directions  /  respectively. 

Therefore,  the  expressions  for  the  induced  velocities  are 
n.  n  p®~**-° 

.  ^my  i*|  fr*tSmC8«+^)-sinCef4-»Vp>-ie^-ep)Q>«(.<^1 


2  $  r 


Je^-W) 

Q*  t^L  1  -  COS  CG>y i  -  8p  •*•  'I’wk-  4>p )  3 

w  An J  C*r-  ir'cosce^-er*>Ju-4v)  +  ^r'ce^-Gp)1- 

The  following  simple  transformation  can  be  used  to  transform 
the  Cartesian  (x,  y,  z)  velocities  into  cylindrical-polar 
(r,  $,  z)  velocities: 

[Up.Cr]*  C  Up.  6*..]  CoS  t8p+-lpp)  +  t^p. Sy  ]  S  iVlCOp  +  typ  ) 

[  Up. <*$,]=  -[U^-CJ  SinCOp+^p)  +■  [UP-£y  Jc-OA  CGp't-^p) 


zrl  ar*co»ce4,-epf  4w-^p>  +  *V*ce^-©p)  J 


The  cylindrical -polar  (r,  $,  z)  velocity  components  thus  can 
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be  written  as 


3*,*  *00 


[uU]=  5~  Ik  fcrxLs.ViCO~-0p  +  tJ;„-q>)  -  ie^-Q*)C 
7^7 +n  ,  C^-  ar'cos<©;,-ep*¥m-%) -*  *»>-*<. 


cos  Ce^-Sp  +  ^-^p)]  ._, 

. - - - ;  i.  r  .  /1ft 


[<Vej>]*^rIk 


0«--« 

,  ew«=  +«o 


—  —  r  1  nn  tk'J 


>V 


tervC»-cov5teZ-  c6u-QP)&.'Ace.;-ep  v^-^p)3  .  . 

0’r- *t'  <^s  t  ©.V,  -  gp  +  4>m-  i^p>  4-  h.*-  r1-  C9„;  -  ©-  )x]  'S/3- 
J6i*  •  -So 

«  +»o 


r  *  [  i  -  Co  s  (  eZ  -  Qp  ■*  4/^,-  vVp  )  ] 


1  C^-i»'vcos(8n;-©p*.^-4'f)+  *'  r 1  ce^ -  Bp >\] */z 

'd  'MK  -  CO 


If,  now,  we  define  xmp  4  (9„,-0p> 
and 

and  transform  the  variable  of  integration,  remembering  that  ep 
is  constant,  there  results 


m  «  I 

n 


CXwp  ~*  Xmp  Cos  CXmo  +  ^wp)  3 

UrV  sr'cos  CXwf+iJ/„p)  +  k^r'xX, r  1  */a 


dx 


mp 


-OO 

pt+oo 


fOl  •  €&]“  — M.  terxC  1  ~  Cot  (  Xwip+vK^p  >  -  Xwip  SinCXwp-t-'I’^p  )  ]  t 

L  £t4TT  &r-1r‘cosCXw,p+^p)+te’-r->tiip  ]**  *”P 


u-oo 
r  +®o 


^TT 


m»  i 


r*£i-  Cosoc.^p  +  iK.p)'] _ 

^■»y i  ar *■  CoS  CxMp  +  ^p)4felrl  X^p  }  3/* 


^Vr»p 


Self  Induction 


Induced  velocity  components  on  any  vortex,  P  ,  result  from  the 
sum  of  the  induced  velocities  due  to  all  rt  vortices.  In  each 
of  the  three  expressions  immediately  above  there  is  a  term  on 
the  right-hand  side  corresponding  to  m.=p  .  This  term  gives 
the  self-induced  velocity  components  for  the  p1K  vortex.  How¬ 
ever,  for  yyi-f  , 

MVp  =*p-4p«o 


and  at  xmp=o£he  denominator  becomes  a  zero  of  higher  order  than 
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the  numerator,  resulting  in  a  singularity  in  the  line  integral. 
The  difficulty  is  fundamental  in  that  the  Biot-Savart  law 
is  invalid  in  a  vortex  core,  i.e.,  where  viscous  forces 
predominate. 


Removal  of  the  Singularity 


If  we  assume  that  the  vortex  has  a  finite  core  of  radius 
the  singularity  in  the  Biot-Savart  integral  can  be  removed  in 
several  possible  ways  developed  to  a  considerable  extent  in  the 
literature.  S.  C.  Crow  [Reference  17]  in  his  work  with  fixed 
wing  trailing  vortices  eliminated  the  singularity  by  cutting 
the  integral  o££  at  some  arc  length  d*on  either  side  of  the 
point  where  lPT„,rl  =  o  .  He  then  subsequently  determined  ct  by 
taking  it  proportional  to  the  radius  £,  of  the  vortex  core  and 
evaluating  the  constant  of  proportionality  by  reference  to  two 
similar  problems  whose  solutions  are  known  by  other  means.  A 
similar  technique  has  been  used  by  Hama  [References  18  and  19] 
for  determining  the  progressive  deformation  of  vortex  filaments 
with  initial  curvature  without  and  with,  respectively,  several 
kinds  of  initial  perturbations.  More  recently,  Widnall  [Refer¬ 
ence  16]  used  the  method  of  matched  asymptotic  expansions  to 
remove  the  singularity  in  the  self-induction  integrals  by  re¬ 
lating  the  problem  for  helices  to  that  for  two-dimensional 
circular  vorte^  rings.  Yet_another  method  is  to  replace  the 
denominator  [I  by  [  |^Pfll-+ 1']**;  this  approach  was  used 

by  Levy  and  Forsdyke  [Reference  15]  in  their  early  work  and 
more  recently  by  P.  C.  Parks  [Reference  2D], 

For  this  study  a  procedure  similar  to  that  of  P.  C.  Parks  was 
selected  because  a  comparison  with  the  Widnall  results  would 
provide  some  insights  into  the  equivalence  of  different  ap¬ 
proaches  and  because  complications  are  avoided  in  this  approach 
when  multiple  interdigitated  helices  are  considered.  Thus, 

C  will  be  replaced  by  e1]1'1  when  self-induced 

velocities  are  calculated  and  subsequently  in  the  perturbation 
equations. 

Even  with  the  treatment  of  singularities  described  above, 
another  problem  remains  in  the  process  of  numerical  integra¬ 
tion,  because  the  integrals  behave  almost  as  though  they  are 
singular  so  far  as  numerical  evaluation  is  concerned.  The 
evaluation  and  interpretation  of  integrals  are  therefore  done 
by  rewriting  the  integrand,  adding  and  subtracting  Junctions 
which  have  the  same  singularities  numerically  but  which  can  be 
more  easily  treated  analytically.  This  procedure  will  be  dis¬ 
cussed  further  in  chapter  5  in  connection  with  integration  of 
the  perturbation  equations. 
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2.4 


VORTICITY  TRANSPORT  THEOREM 


In  this  analysis  we  have  assumed  a  fluid  of  uniform  density  and 
zero  viscosity  and  that  elements  of  the  vortex  line  move  with 
fluid  particles.  If?  denotes  the  position  of  points  on  a 
vortex  line,  the  total  derivative  (Eulerian)  of  ?  is  equal  to 
the  fluid  velocity  at  that  point. 

If  one  is  considering  only  those  fluid  velocities  induced  by 
the  presence  of  vortices,  then,  mathematically,  the  vorticity 

transport  theorem  can  be  expressed  as  ^  C? )  =  IT  for 
points  r  on  the  vortex  filaments. 

A  description  of  vorticity  transport  theorem  can  be  found, 
for  example,  in  the  fluid  dynamics  book  by  G.  K.  Batchelor 
[Reference  21]. 


2.5  PERTURBATION  EQUATIONS 

Let  the  point  on  the  vortex  helix  be  perturbed  to 

the  positions  6x„  ,  y„,+  sy,*  Here  6*^  Sy*,  and  &imare 

considered  to  be  first-order  perturbations;  i.e.,  their  products, 
squares  and  higher  powers  can  be  neglected.  The  functional 
relationships  for  these  perturbations  can  be  written  as 

&*■>»<  =  S  X-rrt  ) 

Sfcwi  «■  fcfrrw  t ) 

where  t  is  the  time  variable. 

Now  writing  =  Cxm  ,  ym  ,  t„) 

and  C6xw,%ym,s.i^) 

the  vorticity  transport  theorem  thus  gives,  for  the  perturbed 
P1"’  vortex ,  +  )  =  UP  +  6tP  > 

where  U^,  (.ip  +  &q>)  is  the  induced  velocity  at  point  rj, 4-srj,  on 
the  perturbed  p"'  vortex. 
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The  induced  velocity,  Up  ,  is  obtained  by  using  the  Biot-Savart 
law  at  the  point  due  to  the  system  of  vortex  helices 

as  modified  by  the  perturbation  displacements. 


U.(.f>+sr.)=  V  -Ca  f  dl:m 

4-tr  wp 

m-|  J  1  ""M 

where  Km*  *  rM  -  ^  +  Sr„^  -  ' 

ciLm  -  d.Lyn  (.r^  v  SrM  ) 

*  [<Ax^,vdC6X«)  ,  dy^+diaa^),aC4-cl(6i^)] 
The  zero**1  and  first-order  cross  products  are 

[R^xcLu]  *  (&-^>x  §-  (?*,'  )de„ 

r  o 

L  Rivip  *  ^  =  (rm  —  Yp)x  &  rm  )  d.©^ 

4-  -  &Vj,  )  X  ,  (.fan  )  d0rvi 

®0i« 

We  may  now  write 

I  Kmp  fb  =  +  gx^-Sx*)1*  Cy^-y„+  &y^-syp)*- 

+  c*«-ip  +  3’3/X 

Defining  the  zeroth-order  term  expression  as 

iRrnpl,,'  -  [CX^-Xp)l+  CvC,-y*)'+  Ci'W-4p)‘L  l”573- 

*  J*PV1 

and  the  first  order  term  expression  as 

I  *  X*lp  C1  *■  ^  Ji*»p  ^  wi  -  Xp7  (.6X.|vt  —  SX.p) 

4  -  yp)  i  syM- ep)  (.64^-  Se^)] ] 


Vi 
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and  expanding,  retaining  only  first- order  terms,  yields 


I  ^rwpl  ,  *  "  Ki*|3 

where  K is  defined  as 

Kmf  =  Cxm  -  xp 'i  C6xm  -  sxp)  +-  ca^-yp)  c  &9™-syP) 

+  (t«  -4p)(6*m -it}) 

The  corresponding  zero  and  first-order  induced  velocity 
vectors (  denoted  by  0  and  1  subscripts)  can  be  written  as 


[upJ#=  ^mP  [<^n“^p)x  )J 

W e|  p  4 1  M 

•m* -°0 

r  e^,»4-«o 

Jmp/4  C<  £i'  - Tp )  X  & ,  aC )  +  (&,  -  Sr,, )  x  i  )] 

,  °®tvi 

;em<-oo 
'  ft*"*00 

K*MpXrM^>  C  ~  Tj,  )x  (fVvi  )J  cl0rvi 


Da 

A-V 


n  „  r 

♦z&-» 

">•>  L. 


It  will  be  convenient  to  carry  out  the  calculations  in  terms 
of  the  helical  coordinates,  r  ande  .  Since  the  parametric 
equation  for  the  m*1  vortex  helix  is 

r Cos  Cemv^) 

y^=  TS(n  ce^  +  4u) 

^r®rC, 


the  first-order  displacement  perturbation  expressions  can  be 
written  as 


&x^  =  srw  Cosce«*- -  r  5m  t©M  +  ^)6(t>w 


&\C  =  6rw,s.n(©ft;+^) +-  rcos  (.0^+^)  &<t>, 
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Sim  *  Sim 

The  zero**1  and  first  order  terms  of  the  numerator  of  the  inte 
grand  in  the  *  ,  y  and  a  directions  also  can  be  written  as 

Zero1*1  order  terms  in  the  x  direction 


■  yp  )  ^L,  (im  )  -  (4m  "  if  )  ]L,  Cam  )  *  f  T  Sin  (Sm  +  'J'm)  -  T  Sin  (0^+  )}  •  kr 

—  hr  (,6m  "  0p)  CoA  (  Gm  +■  'Kyi ) 

First  order  terms  in  the  <  direction 

-  (Sim -Sip)  3  = 

^0m  36m 

Rr[5fm  SiY)(0m+^)  -  SmCepV^)  4-  CoS  (6m  +  'I'm  )  -  fS4,pcosCV'H0] 


-  CSim-Sipl.rCosCeM  +  ^m) 

(ym-yP)  ^,  (6-iw)  -  C6tfrn)  - 

*6 m  66m 

r  3_J  (Sim)[Sin(e^  4-'KWi)  -  S»w(.6p4-IJV)J-  (6m-©p)Cfs'^6r'4i^r«(©i»1+vVm) 

*■  6r„Cos(.Q^*-^)  +  rCos(0Mr^]  ®  ,  C6c>m)  -  rS^SinCe^  +  i^v,)] 


Zero**1  order  terms  in  the  y  direction 
Cim-  ip)  3  (Xm  )-  (Xm-Xp)  ,  (4m>  * 

-  hrtcef;-ep)SmC0^<-4'w,)  -  kr^CosCdm  +  ^m)-  Cc&C©p+-'Vp>} 
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F irst - order  terms  in  the  v  direction 


dim-  6*,)  a,  -  C6X*-*Xp)  (.*'„)  - 

SiV»(e„  +  il»m)  -  ftrCtr^CosC©^*^-  Srp  Cos  C©p+-4*p ) 
-rs^SiflCe^ +  *64>p  Sm  (ep-*->VP)] 

Rr(e^-©p)  [|lb<  isr-m)  cos(en;+4'w,)  -  Sr*,  Siwce^+v^)- r6$fV,Cosle^+'U 
-rSin (e>4u)^,C6<t>J  “  v&.CMwOfrww^MW-  CosiGp+^p)] 

0*TV1 

Zerot^- order  terms  in  the  a  direction 
cxi-xp)  -  cyj,-yp)^cx^,)  - 

r'cbs  (erl,+'ky,){Cos(«^+4»*,)-  Cos(©p-v^f)j.  _  r's  Sm(p^} 

First-order  terms  in  the  »  direction 
c6jC-8Xp)^ca^)-(sy^-syp)  ^t*m>  - 

rCoS(€C  +  v|;m>[6rmCos(ew+^m)-  Sfp  Cos  £©,.■»■  ^p>  -  r&^^Sivi  0^,4- l^ry,) 

+•  YS^pS«n(ep*-ll/p)j  +•  rS»r»C©m  +  »(\wl  [.Sr*,  SinC©^*-^) 

.  srp  s.v.cep+  \\JP)  +  rsfc,*,tosc©n>4>wO  -r^s.vtCepv^p)  3 
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r|cascam  +  4im)  -  c©*cep*4v>]-L&;,(6r'">  Sir»(firA+4u>+  6»v* Co* ce^*4»»0 
-rsdi*  rtos(e^+4'rT.^/C6<t)J-'r|.Since^+^-s.v>ce^4.^)] 

[  ij;srM)Co*(e^*-ilu.i  -  Sf*  Sin ts^+'Vrwi)  - TS4>wcosi©^+4-w)-rS,n(e^^^6^ 

Jm\*  ar- av‘*-Co*Ce^-©p-»-4'w,-^>  -*-fclr'te^-ep)’' 

K*ip  -  rSrM{i-coi(x^4-^f)}  +  r6rp  {  t-CosCx.mp+^p>] 

♦  ^  Sin  CXw.p4-4w^>  -  r'6$p  &mCx^p4- 

^  Xn«p  6tro  ter  Xrvtp  6^^ 

Calculation  of  Terms  on  the  Left  Side  of  the  Vortex  Transport 
Perturbation  Equation 

sxp  s  r  Cos<e^-*-4»p)  +  srp  costo^-t  q>p)  -  rs»r*cep-Mj>>  &$p 

• 

^cxp+6xp)  r  rpCoscep+(4»p)-  rfy  *  sfpCos  (ep+vj,p) 

-  6rp,  4>p  S,Vi(©p+4>)  _  fp&4>p  Sm(6p  +  lVp) 

■  • 

-  y4>P&<tJpCosCep*-'4't-o  -  rs<J>p  SinOpf^)  (2) 
a^+  6yp  =  r smOp  +  v^p)  -*-  s.vycep+^fO  +  y&4>pco&f©p  +  ^p.) 
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icap  +  sy*)  =  Yp  S»nCep  +  4>p)  +  Cos  i©p  v  )  +  &fp  S>vt  (.©p.H+p) 

+  &Tp4>^CoS(ep+4'p)  4  ^p&^p  ^OS(0p  +  4/p) 

+  r&^pCoscep-MVp)  -  r^,  &4>p  s>v<  cep+ 4>p  >  (3) 

^  +  &*>  (4) 

We  now  resolve  the  velocity  radially,  circumferentially,  and 
axially  (i.e.,  in  r  ,  <£  and  4-.  directions) ;  these  components 
will  be  defined  as  i  ,  rf  and  ^respectively. 

K  p  r 

Multiplying  equation  (2)  by  CosCvV^A  and  equation  (3)  by 
Sm(VV^)  and  remembering  that 

Cos<op+^p+  )  a  Costep+q/p)  -  &4>pSmcep*-4/p) 

SiV,  CSp+ l|>+  6(j>p)  ~  StV,  CSp  +  ^p)+  6^pO>^(0p+ q>p  ) 

we  obtain  the  total  radial  velocity  of  the  perturbed  helix  at 
p;  thus: 

^pc  +P  +  &*>“  T4>p  6<J>p  +  TCjjpfc^p  =  TTp-PfeVp 

where  f>  ■  zeroth-order  radial  velocity  at  P  . 

sifp  ■  first-order  radial  velocity  at  P . 

Similarly,  multiplying  equation  (2)  by- Sinfe^^+fcijjp)  and 
equation  (3)  by  CosCep+qp+tfy)*  we  obtain 

f[p  =  r4>p  -  rp  s4>p  +  srp  cj>p  4-  rp  s4>p  -*■  r  &<j>p 

=  +  &r^+  r  S*^ 


18 


*  aL 

where  f<ppis  the  zero  -order  velocity  and  the  remaining  terms 

on  the  right-hand  side  represent  first-order  circumferential 
velocity. 

Finally,  ^  -  ip+**p 

where  3p  *  zero^-order  axial  velocity. 

fcip  *  first-order  axial  velocity. 

Transformation  of  the  Induced  Velocities  to  Polar  Form 


The  zero^-order  velocity  in  the  x  ,  y  and  t  directions 
which  result  when  the  right  side  «~r  equation  (1)  is  expanded 
will  be  called  U*p,  Uyp  and  Dtp  ,  respectively,  and  the 

corresponding  first-order  velocity  terms  will  be  called  fcu*p,S>U3f 

and  &u*.  Expressing  the  induced  velocities  also  in  terms  of 
radial,  circumferential,  and  axial  components,  and  equating 
the  left  sides  and  right  sides  of  equation  (1)  component  by 
component  yields 


ip*  +  =(Uxp  +  &l)jl.p)  Cos  (6p+  p ) 

SU*p)  StVi  cep 4-  4*p+  «4>p) 


%*r<pp+  Srp4>p  +•  r 6<J)p  *  -CoKf+SUxp)  S(n  C©pf%+  Sffcp-) 

-t-  c^p  +  sUap)  CosCep-*-4’p+'6<t>FO 


Now,  separating  zero**1 -order  and  first-order  terms,  the  vortex 
transport  equations  are  obtained  in  terms  of  helical  coordi¬ 
nates  for  unperturbed  and  first-order  perturbed  motion,  as 
follows: 


Zero^  -order  (unperturbed)  transport  relations 
"fp  *  U*,p  Cos  (6p+l)>p)  +  Uyp  S»V\  C8p+  4*P  ) 
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t^p*  -  0»p  SirCep+'Vf)  Ujf  Cos£ep*4v) 
ip-  U*p 


First-order  perturbation  transport  equations 

“  Uxp  S$p  +  +■  Uyp  i(^p  Coi<©pf  4^p) 


+•  SUxp  Castep+^P)  +■  SUjjp  Sin  C©P  *-4»p) 

*  T<j>p  6$p  4-  &0xpCosC©p4  typ)  4-  SUapSiYlCSp+^p) 
6rp<j>f4-rStyp  =  -  U*p  %4>pC<M(ep4^p)  -  Ujp  &$p  Simcepviyp) 

-  S0*p  SinCOpf +p)  *  &Uyp  C&*  Cepv  4*p) 

*  -  Yp  t$p  -  &Uxp  Sim  C6p4-  p  )  +  SUyp  CosC©p4'i4>p) 


or 

srp  =  r4>p^r+6Uxp  cosiepv  o>p)  +■  6i)y p  Since p*^p)  (5) 

-%rr4y-fp6«^f  -  SOxf  Smcef+«|»p)  +  SUyp  Co*  C©r4lW  (6) 

sip  »  S0*p  (7) 

Substituting  for  induced  velocities  and  time  rates  of  displace¬ 
ment,  the  unperturbed  relations  become 


r,.y  * 

r  Z_  4ir 


hr  ^ Sin  (Xmp 4-  -  y rwp  CoS  (X***p  4- tymp)3  Jiwp  dXmJ> 


n  r°°  ^  3 

4>p« ^  Oil  I  i CoS  (XiMp4  *Vimp)  -  Xr*if>  Sim  CX)vip4-4'mp>  J  Jmp  dLXmp 
-to 


n 


ATT 


me  I 


T  ^  >  -Cos  (H*vip4-'4'mp)  j  TVwp  4jCrvip 


-00 


20 


Similarly,  the  first ‘order  perturbation  equations  can  be 
written  as 


-  kT  Xmp  CoS  (Xmp4  V  kr  Sfrn  Sitl  (>C**<p4-  'J'Vyip)^ 

+  *■  Xmp&4>MSlnUw,pv4>mp) 


“  X,,,P  A/ifidim)  Co*lx»p+4»iwp)  -  8<bpCos(Xwp+ 


+  r  ^-SimCos(XwpnV,*p)  *  &tpCo»(Xrwp4-lVmlO 

+  J|^<S**0  Sin  Cx»p+-  4’i*ip)J]clxmp 
A 

~*7~  ^pp1  hr'^SmtXmp  +  i|/»ip  )  -  XmpCostXmp+^p)}  kmf  dxMp 

"•»  *  J 

\  (8) 

«  r*° 

■S^»  ^^.Ga  I  Tmp*‘L^rl,t^p^^6r»*i^G»»CXiMp+-4»iwp)-6fmG9Stxw,p+-vVMtp) 
**i»i  J 

~  X*wp  STWi  SiV>  CXmp-t-  l^mp)  +  6Yp  CollXmp+-  l|Vnp) 

+  Xmp  Bfp  Sin  lXmp4»  l|>iv»p)  J  +  hr*'^-X,*ip;ft.((6d^|>Sin0lmp*-((tip) 

~XiMp  %4>mCo*(Xmp>.  U»mp) 

+  *mp  Sdp  CoSlXmp*  ^mp)}  +  lr^,"  ^e»vi  SmCX»*»p+  V|>mp> 

+  Sip  &tV\ (. Xtwp +■ l^mp ) 


n  r° 

^Jr  kflJWp  *.  Kr*ip  [i 
mti  J 


1-CosUmpt  ^P  - X^p Sin  tx  *,pf vV^p^cIjChp 
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felrp  a  \  *  r* 


^_4ir 

W«  I 

-co 


«U 


Jmp1  [r  {|^(C6n«)  s,vx  oc-pvt^p )  3.  hr, 


M 


—  Cos  Cx-^p-v-l^vip  ^  ~  &Tp  Cus  (X-ivip^"  4<*ip)3 
*  rx{jL<64,)  +  S^SmCx^p+^-Cos^^^/S^ 

-^PS  ivt  tX.*yip  +  Vwp  )  3  j  dx-kwp 


n  r® 

~^~  Og.)J^,r/it'r'1{ i -cos cx^p+.iiUp)}  ^Kv»p  dx.* 
Mei  J  on 


CIO) 


=  t  -i-w 

Removal  of  Singularities 


As. explained  in  Section  2.3.2,  the  integrand  that  provides  the 
self- induced  contribution  to  induced  velocity  becomes  singular 
at  Xmp  ■  0.  This  singularity  is  eliminated  by  changing  the 

expression  for  JV^p  from  [ ir- zrxCot> C.Xwp+  4iwp)+  /,_ 

to  ijr’c.o*<:>tr^p  +  ^nflf )  +  -*•  ^  •  } 

when  **>•■  P  .  Here  £.0is  the  vortex  core  radius  and  is  small 
compared  to  the  radius  r  of  the  helix.  Note,  however,  that 
when  wij<p,the  expression  for  Jmp  remains  unmodified. 


2.6  TYPE  OF  PERTURBATIONS  AND  EIGENVALUE  EQUATIONS 


The  integro-differential  perturbation  equations  derived  in 
Section  2.5  adi^it  solutions  of  the  exponential  type,  for 
perturbations  sr^  ;  that  is, 
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eKt:t  t"9"' 


(U) 


or  in  component  form 


S>r*.  = 


tfyv,  €. 


*t+-  ttO©^ 


.  ^  «( t  +  L  lO 

Sk*  s*„e 


*  ort  v  cujem 

S2*= 


-1h 

where  Sr„,  -  vector  perturbation  of  the  vn.  vortex  having  sfM, 
S^and  S^^as  its  components  in  the  r  ,<£  and -a- 
directions . 

oC  ■  exponential  growth  rate  factor,  complex. 

co  ■  wave  number*of  the  perturbation. 

In  principle,  no  generality  is  lost  in  considering  such  solu¬ 
tions,  since  an  arbitrary  perturbation  can  be  synthesized  from 
these  by  Fourier  Integration.  In  practice,  the  numerical  work 
will,  of  course,  deal  with  a  limited  number  of  wavelengths. 

The  exponential  form  assumed  in  the  solutions  of  perturbation 
equations  (8),  (9),  and  (10)  allows  this  set  of  integro-differ- 
ential  equations  to  be  written  as  a  set  of  linear  algebraic  A 
equations  in  the  amplitudes  of  the  perturbation  components  6rm , 

The  spatial  derivatives  of  the  solution  forms  are 

given  by 

i,tsr„)=  1  [&r^eo<t^<o0,v'  1  =  infirm 


*  The  inverse  of  the  wavelength,  i.e.,  number  of  waves  per 
helix  turn. 
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and 


Also,  the  £irst-order  displacement  perturbation  components 
can  be  rewritten  as 

A  i-COSy^  .A  «t  +  LWO^-itoeP  +  icOOf, 

b'nt  «  orWi  «  a  %rm  e 

»  6rw,e‘k  +  iu,®p 


6<J^  =  c«t+  twQp  gitoXmp 


tew  el~x"P 


(12) 


Substituting  equations  (12)  into  the  perturbation  equations 
(8),  (9)  and  (10)  and  defining  4.  Xmf  •«■  chnp  yields  the 
following  eigenvalue  equations: 


4TT«  srP 
p 


hr6fw,  Sir»y*irtp  Sin  toxmp 


-  Xrnp  Co*  9mp  CosCOXmp  -V  Sin Jimp  Cos  WXmp }  4-  i  fcr&fm 
(-to  Xmp  Sin  CoS  coXrvip  —  X^p  CoS  bimp  SinCOXmp 

+  Sin  Simp  fcincOXmpJ  +  ftr1'6^w>^Co*ympCo*W>Xmp 
+  Xmp  Sin  CoS  CO  Xmp  +  CO  Xmp  Cos  y*ip  S>n«OXmpJ 
+•  iler^fc^^^Ceaymp*'^  c°><",iy'P  4-X.ynp  S'V>9mpi*vsOM«p 

-  coXmpCosy^pc««  eoxmpj  +  rS*»v.t  -  CosSJmpCoS  «*«p 

-  to  Sin  Sift  CO  Xmp J  +■  lv&%  m  ^  -  CoS  ^mp  S'  m  MKwf 

%  ^  X 

+  U3 S'Vi y m p  Co-S  J  +-  XV  6^PV’CoS^,**f-^mpSin'JiwfJ 
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r  Sip  cos  Xw,p^  dx.,vip 

Jj,^pS/t[WrV{ S'V'y^p  "  *-,v1p  006  ^p}  J .  [r 

-"{(,-CoS  )  CoS  CO  -4-i(l  -  Co4  y^p)  Sl -JO  to  M-^p  jr 

+-  r*s4>w  ^  SiViy^pCas  cox-^p  +  1  s.v*  'd^ps'v'  «*~p  } 

4-  ftr  Stw  {  4  x«,  coi  (OX.A,p  +  L>o~t  wx-mp] 

+  v-  6rp  ^  i  -  cos  y*v,p]  s>  ^  '^y'/'f>l  “  ^ 


Similarly , 


(13) 


OO 

^MTf°t  S»P  m  fj^p1/x[ter6r^  ^-cox.wpCosy^pS'V'usx^p 


rV»»|  -°0 


-  Cos  y,v,p  CoS  cO^'V’p  —  Xyv>  p  SiVl  Hivip  CoS  OOXwip^ 


+  Ifer  6  rw,  ^COXwip  CoS  y  rn  p  Cos  LOXwnp  —  Cos  y^p  S\YKoX*v,p 

-  JCwpSmyWpSiV  ooxrvip3  +  Rrfirp  ^Cjai/rv,^-x«'p,'y'^’*'p} 

4  terlS4jv,\u3>Uv1P  SiV\y**p  SiV>  OOXn^p  -X-nnpCosy^p  Ca&cOX-wp 
4  5iY»ywpCoS£OX»vip]  +il«>'l6^M{-w,CwP  SiVjy^pCo-SUJK^p 

-  Xw,p  Casy^p  SmcOXwip  +  SiVi  p  SvVi  WX«,p  3 

-  ^sjptsmy^p-  Xw,pCosy^p}+  rS^-Sfvi^^^p 
4-  coCosy^pSfvi  u>Xv*p  -  co  Si-n  ooX-w'p  3  4-  t-rS  ^ 

{-  Sm^p  cOXwip- cocaSy^pCostOXwip  +  coCoSCOXw,^ 

■+■  rfitp  Si'ta  y*^  p  0  dxwip 


-  3 


a  r  c 

y  j^p  /z[*r*{  I  -  Cos  ymp  -  Xmp  S.V« 

wir  I  J-«0 

{  4- (I- Cosy^ip)  Cos  COX»y,p  4  C  Cv-CosWwip 


}]■  LrS^ 

)  S\v\  LOX^vip  ^  +"y  64>m 


{  4  Si'ioywpCa6  COXvvtp  +*  ^SiviympS|,'VVJOX.m  p]  4-  RrSfcyvi 

-  'V 

{4  Xrnp  Cos  CO  Xy«p  4  L  Xyv,p  SiVi  OJXyvip  *J  4-»'S»p 

{.4-  1-Cosy**^  4V  S4j,^S"'><!/*vip3 -fer'6tpX|vip  J  ci  Xmp  - 
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^  CO  S  iv')  y*v»p  StV>  UlK.^ p  4-  2_  CoS  COXwn  p  —  Cds  cj^  pCos  tOX^pj 

-CO 

+  CoSt^^y^pCoSuDK^p  4  ^S^y^tO>L^p  — Cosy^pS-ontO^p?) 

A  £  r 

—  Y  Srp  Ccxsy»v»p  “V  Y*  S  v-  c-Oc5 m  l>0  X*v*p  -4-  St'-o  v^3  kna  p  CT<j»s  u)X  yv^  p 

A 

4-  CO  Cos  p  S  *'v\  LO  X.  W*i  p  5  +  \  CO  Cos  CO  X-v-v^p 


4-rrocfe  ip 


V\ 

=y 


-+■  SiV\  y ^  p  S iVi  co  )Cw p  —  co  c©a  si »vi p  Ccs  co^wp  } 


-i/a. 


■+■  y  ^  S  ^  ^ ^  4  w\  p  ]  ]  a  Xw,p  •  J*rvy  p 

r°° 

"3Z  >].L  Y&r™  l  C-cosy^ 


Cos  £OX**p 


-oo 


4-  1  (.  '-  c^»-s  y^p)  S\Y\  COKv«p  ^  ■+  V'l&cL^ 

■+  St'vi p Gci co +■  t  £v'v>y*v,p  5>V>  w>iMp  ^ 

A  C  A 

4-  ^y^i^^XrvAp  Cos  CO Xw> p  -4-  c  M~vv»  p  Sv‘vy  COXvn  p  j  4-  V  S/p 


{'-cosiinp}  -  r1- 

*  1^  r  Sip  X.w\p  cl  X-  p 

2.7  STABILITY  ANALYSIS 


Equations  (13)  through  (IS)  can  be  written  for  integer  values 
of  p  ranging  from  1  to  in  ,  where  re  corresponds  to  the  number  of 
interdigitated  helices.  After  numerically  evaluating  the  set 
of  self  and  mutual  induction  integrals,  a  system  of  3 n.  eigen¬ 
value  equations  is  obtained.  The  limits  of  integration  are 
from  -<o  to+to.  The  integrals  in  the  perturbation  equations 

-^2.  ~  S/i.  -3 

are  convergent,  since  both  J^p  and  Jrvip  vary  as  X^p  and  X,y,p  , 
respectively,  in  the  limit  as  approaches  infinity. 
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Nevertheless,  these  integrals  must  be  calculated  numerically 
with  some  care  and  with  proper  treatment  of  the  singularities. 
These  techniques  are  discussed  in  the  next  chapter  in  more 
detail. 


2.8  SYMMETRY  CONSIDERATIONS  FOP. 


For  two-bladed  rotors,  it  is  possible  to  separate  symmetric 
and  antisymmetric  perturbation  modes  since  the  resulting 
eigenmatrix  is  symmetric.  No  such  simple  reduction  in  the 
order  of  the  problem  can  be  made  for  rotors  with  more  than 
two  blades,  as  the  symmetry  of  the  eigenmatrix  is  lost.  The 
procedure  is  as  follows:  define  *,  and  as  the  perturbation 
vectors  for  the  first  and  second  helix,  respectively; 


that  is, 


'  *  1 
Sr, 

r 

Srv 

A 

,  r6t, 

and 

tx 

r» 

H> 

A 

/V 

-  -> 

- 

and  [x,]85  the  perturbation  vector  of  the  matrix 

equation  Ax  Ax.  Equations  13,  14,  and  15  show 

that  At 

The  eigenmatrix  A  can  be  written  in  the  form 


where  A„  and  are  3x3  matrices  and  represent  self  induction 
for  the  first  and  second  helix,  respectively.  and  are  the 
3x3  matrices  representing  the  mutual-inductance  effects. 


The  eigenvalue  equations  can  be  expanded  in  the  form 


>  x ,  =  Anx,  +•  A 12.  Xx 

AXl  =  A 


Consider  a  class  of  perturbations  in  which  disturbance  motions 
in  one  helix  are  repeated  on  another,  but  not  necessarily  at 
the  same  time.  This  kind  of  disturbance  can  be  investigated 
by  imposing  the  relationship  =  e'*Xi.  In  these  circum¬ 
stances,  the  above  two  equations  remain  compatible,  since 
A. i  ■  £  xxand  An.  *  A*»  •  Three  eigenvalues  can  then  be  solved 

from 

Ax»=  e  ^An  Xi.  +  —  CV*dw  +■ 

These  eigenvalues  will  give  maximum  and  minimum  amplification 
rates  for  certain  values  of  ^  .  In  Figure  2,  divergence  rates 
are  plotted  against  the  phase  difference  .  For  |  *  0  or  2n , 
Xi  ■  x*  ;  for  $  -  ir,  X,  --x».  Thus,  these  are  cases  of 

special  symmetry:  in  the  first  case, perturbations  in  the  two 
helices  are  of  same  magnitude  and  occur  in  same  direction 
simultaneously;  in  the  second  case, they  are  of  the  same  magni¬ 
tude  but  occur  in  opposite  directions  at  the  same  time. 

These  special  cases  have  been  incorporated  in  the  computer  pro¬ 
gram  for  the  increased  efficiency  afforded  by  that  option.  The 
full  eigenmatrix, which  does  not  involve  the  special  symmetry 
combinations  described  above,  has  also  been  calculated  for  the 
two  helix  cases,  and  the  results  are  identical. 
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CHAPTER  3 


COMPUTER  PROGRAM 


3.1  GENERAL  DESCRIPTION 


The  computer  program  for  the  stability  analysis  consists  of  a 
main  program  and  integration  and  eigenvalue-analysis  subrou¬ 
tines.  The  main  program  has  three  significant  subparts,  which 
can  be  described  as  follows.  The  first  part  deals  with  treat¬ 
ment  of  singularities  for  the  self-induction  integrals.  The 
second  part  deals  with  integration  and  with  obtaining  the  per¬ 
turbation  equations.  The  third  part  is  that  used  for  eigen¬ 
value  and  eigenfunction  determination  and  printout.  The  three 
subparts  of  the  main  program  and  the  subroutines  will  be  des¬ 
cribed  separately  in  more  detail  in  subsequent  sections. 

The  program  is  written  in  FORTRAN  IV  language  and  has  been  run 
on  an  IBM360/65  computer.  Double  precision,  i.e.,  operation 
with  16  significant  digits,  has  been  used  in  the  integration 
and  in  the  portion  of  the  program  eliminating  the  singularity. 
The  rest  of  the  program  and  subroutines  are  in  single  preci¬ 
sion,  i.e.,  operations  in  which  8  significant  digits  are  car¬ 
ried.  Wherever  a  "substitute  function"  approximation  has  been 
made,  numerical  checks  have  been  made  to  insure  that  a  more 
complete  approximation  would  not  change  the  eigenvalues  to 
within  the  last  5  significant  digits.  A  block  diagram  for  the 
computer  program  is  shown  in  Figure  3. 


3.2  TREATMENT  OF  SINGULARITIES  IN  THE  SELF-INDUCTION 

- rmrmrs - 


As  explained  in  Section  2.3.2,  the  singularities  in  self-induc¬ 
tion  integrals  have  been  eliminated  by  considering  that  the 
vortex  cores  have  a  finite  size.  That  is,  the  term 

Jr/1*  [^-«,c«(Kfp+  ]'3/t 

has  been  changed  to  TPpV‘ =  [ar v-  ir,to*«^+  vy^) +  hV-xPF>  +  . 

Although  classical  singular  behavior  of  the  integral  has  been 
eliminated  by  this  modification,  the  integral  is  still  quite 
sensitive  to  numerical  evaluation.  Figures  4  and  5  show  some 
of  the  integrands  which  occur  in  the  perturbation  equations. 

The  former  exhibits  near  singular  behavior  at 
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As  noted  earlier,  numerical  evaluation  of  the  integrals  has 
been  accomplished  by  adding  or  subtracting  functions  which  can 
eliminate  the  troublesome  near-singular  behavior  of  the  inte¬ 
grands  and  which  can  themselves  be  integrated  in  closed  form. 
The  following  are  eight  such  integrals  which  have  been  found 
useful  in  dealing  with  numerical  difficulties,  together  with 
their  analytically  integrated  values. 


Denote 


a.'1  i 


0  +  *1) 


and 


*  £  Sinh'(iy) 

Q-  ' 


where  xu  *  some  upper  limit  to  be  chosen. 


I- 

XU 

T  =  r  xdx  _  J _  fj _ _ ! -  1 

^  Jc*AcOOv>-~  Ci4-y»'lLt 

O 

*U 

o 

fXU 

I  *  U-HO'^U-X:  =  -J - -  [tcuwlAt-  ifccuv^i] 

C 

I_«  (i+  hx)'5/»-(  )CcLk  _  _  J! _ -  f—  -  - - - - — v  3 

JlX,-+av>5>,t  3ct+-fe\)  L  E*  \ 

o 

rK0 

I,  =  =  J _  [ tcxnh3*! 

4  Jck1+<*.x)s/‘-  ja1 

o 

r*u 

T  r  (  x>ot* -  =  — _  r^_  — _ _ _ _  +■  - - - -] 

rxu 

X«=  '*■  *■*  ^ x - -  ~J -  -  -tanh  -  i  tcxviW* -t  1 

■  0+-W\)x 
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In  Figure  6,  the  dashed  curve  II  shows  the  plot  of  one  of  the 
singular  integrals  plotted  after  subtracting  these  analytical 
functions  in  the  range  0  tc  stt  .  Since  *.0  is  an  arbitrary 
upper  limit  on  the  above  integrals,  it  is  chosen  such  that  the 
addition  and  subtraction  mentioned  above  will  be  carried  out 
only  in  the  region  w'here  the  original  integrand  is  nearly  sin¬ 
gular. 

It  can  be  seen  that,  for  small  values  of  xpp ,  the  modified 
function  ~  [iy1- ir lo- s-rt )  +•  x-pp  v  to  3 

2  C  c  r'  VEc1-]  =  CH-tev>L  r'-x-pp+  3 

Similarly,  if  in  the  numerator  one  makes  the  apnroximation 
cosXppSi  and  SiViXpp=  xpp,  the  perturbatirn  integrals 
take  the  form  of  one  of  the  integrals  I,  to  .  Therefore,  the 
perturbation  integrals  behave  like  integrals  I,  to  I8  in  the 
neighborhood  of  *ppB  0.  Some  idea  of  the  errors  involved  in 
this  approximation  is  gained  by  the  comparison  of  perturbation 
integrand  and  modified  perturbation  integrand  plotted  in 
Figure  6. 

The  analytical  integrals  used  to  remove  near-singular  behavior 
labeled  I,  to  I*  and  mentioned  in  the  previous  section  have  been 
denoted  by  symbols  coco  to  coc8)  in  the  computer  program.  The 
integrand  of  the  perturbation  integrals  is  thus  modified  by 
subtracting  coco  to  coca)  from  the  appropriate  integrands.  The 
upper  limit,  xu  ,  on  the  integrals  T,  tor8  has  been  chosen  to 
be  2tt;  beyond  this,  the  perturbation  integrals  behave  suffi¬ 
ciently  smoothly  to  be  amenable  to  standard  numerical  tech¬ 
niques. 

M  is  the  parameter  which  selects  the  overall  integration 
limits  on  the  perturbation  integrals;  M  ■  1  corresponds  to 
semi- inf inite  integration  limits,  while  M  =  2  corresponds  to 
doubly  infinite  integration  limits. 

The  process  of  integrating  the  modified  nonsingular  expressions 
was  tested  by  calculating  some  sample  integrands  against 
increasing  ,  to  determine  a  reasonable  cutoff  upper  limit 
for  infinite  integrals.  This  limit  is  dependent  upon  the  pitch 
of  the  helix.  For  a  pitch  of  less  than  0.2  radian,  the  inte¬ 
gration  has  been  performed  up  to  IOOtt  =  314  radians.  For  pitch 
values  between  0.2  and  0.4  radian,  the  cutoff  upper  limit  was 
set  at  70tt.  Finally,  the  cutoff  limit  for  values  of  pitch 
ranging  between  0.4  and  1.0  is  50ir.  Increasing  the  cutoff 
limit  beyond  the  values  specified  above  does  not  change  eigen¬ 
values  to  within  five  digits. 

A  combination  of  two  numerical  integration  methods  has  been 
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used  to  balance  accuracy  and  computer  time.  Since  the  major 
contributions  to  the  integral  are  made  when  is  small,  a 
double  precision  Simpson's  rule  [see,  for  example,  Hildebrand, 
Reference  22]  has  been  used  with  an  interval  size  limit  of  14tr. 
Between  14tr  and  the  appropriate  cutoff  upper  limit,  a  double 
precision  quadrature  integration  has  been  used.  The  quadrature 
integration  technique  is  much  faster,  although  less  accurate. 
Parameter  Ll  of  the  main  program  preselects  between  the  two 
separate  sets  of  integration  subroutines.  Also,  the  parameter 
iLi  of  the  main  program  stops  modifying  the  singular  integrals 
beyond  an  upper  limit  of  xu  ■  2ir. 


A  short  note  regarding  the  quadrature  integration  technique  is 
perhaps  of  some  value.  To  compute 


y  = 


fb* 

■fc*.)  cLx. 


Ao 


by  Gaussian  quadrature  formula,  an  a.  point  expansion  in  terms 
of  Legendre  polynomials  is  used.  First  transforming  the  range 
x  ■  a,  to  b,  into  T  ■  -i  to  +  i  by  defining 


ax-ccu-t-b.) 

Cb,-a.) 

there  results 


or 


X  c 


b»-a»  ^ 
2. 


b,-»  a. 
2 


b.-a. 


.+ 1 


4>CO<At 


with 


e  -f  bjrjJU^ 


-I 


Using 


.  in)  ,C«> 

for  some  given  n.  with  coefficients  A*  and  nodes  (note  that 
t*  are  the  roots  of  Legendre  polynomials  of  degree  n  ),  the 
result  is  the  approximation 

yt,*  Cb.-a.)^T  [0*-a.)^"  +  *°±9"]} 

11*1 

which  is  exact  whenever  f(x)  is  a  polynomial  up  to  the  degree 

ln-t. 
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In  the  present  program  the  subroutine  Q&9  evaluates  the  inte¬ 
grals  by  means  of  a  12 -point  Gaussian  quadrature  formula  which 
integrates  polynomials  up  to  degree  23  exactly.  This  subrou¬ 
tine  is  based  on  DQG12  taken  from  IBM  publication  System/360 
Scientific  Subroutine  Package,  Version  III,  Programmer’s  Man¬ 


ual. 


The  nodes 


i.CrV) 

and  coefficients 


used  in  QG9  are  listed 


in  Krylov  [Reference  23],  IBM- supplied  Scientific  Subroutine 
Package  DQG12  has  been  modified  to  collaborate  with  subroutine 
INTG,  which  supplies  the  integrand  functions. 

Test  programs  were  run  to  check  the  convergence  of  integrated 
values.  The  upper  limit  on  Simpson's  rule, integration,  the 
range  of  integration  in  one  step  of  the  Gaussian  quadrature 
formula,  the  interval  size  in  Simpson's  rule,  and  the  upper  cut¬ 
off  limits  for  the  infinite  integrals  were  all  varied  till  the 
eigenvalues  converged  to  five-digit  accuracy  for  the  range  of 
parameters  in  the  program. 


3.3  FORMULATION  OF  EIGENVALUE  EQUATIONS 


The  perturbation  equations  of  Section  2.6  are  a  set  of  eigen¬ 
value  equations  consisting  of  the  growth  rate  terms  on  some 

p*h  vortex.  These  growth  rate  terms  are  dependent  upon  the 
self- induced  velocities  due  to  perturbations  in  the  position 

of  the  p1* vortex  and  the  velocities  induced  by  the  remaining 
(  n  -  1)  vortices.  Note  that  mutual- inductance  terms  influence 

the  vortex  in  two  ways:  first,  by  affecting  the  growth  rate 
of  perturbations  on  the  p‘1hvortex  due  to  velocity  changes  in¬ 
duced  by  deformations  of  the  vortices  other  than  the  p^vortex; 

and  second,  by  inducing  zero^-order  velocities  on  the  p^vor- 
tex  which  affects  growth  rate  terms  involving  perturbations  of 

the  p^vortex  itself. 

The  terms  associated  with  the  vn  vortices  and  contributing  to 

the  growth  rate  of  the  p1wvortex  differ  from  each  other  through 
the  constant  ;  values  of  this  constant  are 
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tabulated  for  each  p  and  m  in  Table  I.  Each  value  of  P  and 
m corresponds  to  contributions  of  the  perturbations  of  the 

vortex  to  the  growth  rate  of  the  vortex.  This  contribu¬ 
tion  is  accounted  for  in  terms  of  a  3x3  matrix.  It  will  not, 
however,  be  necessary  to  obtain,  numerically,  all  nxm  sets  of 
3x3  matrices. 

Since  the  effect  of  some  m"' vortex  on  the  p*' vortex  is  depen¬ 
dent  only  upon  the  difference  rather  than  either 

H^or  4^,  the  resulting  symmetry  considerations  reduced  the 

numerical  work  considerably.  All  the  diagonal  elements  in 
Table  I  correspond  to  0;  therefore,  there  is  one 

such  3x3  matrix,  representing  self-induction  effects  in  the 
three  coordinate  directions.  All  the  super-diagonal  elements 
correspond  to  the  mutual  inductance  terms  [p*l,  rr»  *2 ;  p  *2  , 
^■3;  .  .  p"  n-1,  tn«n]  and  have  .  These  also 

need  only  be  evaluated  once  instead  of  (n-1)  times.  It  follows 
that  in  order  to  obtain  the  complete  system  of  eigenvalue 
equations,  only  (  2n-l)  set  of  separate  calculations  need  be 
made. 

All  of  these  reductions  in  the  number  of  additional  evaluations 
of  eigenmatrix  elements  occur  whether  semi- infinite  or  doubly 
infinite  vortices  are  considered.  Further  reductions  in  the 
number  of  calculations  which  must  be  performed  are  possible  for 
the  doubly  inf inite  case,  i.e. ,  whereM-  2,  and  such  simplifi¬ 
cations  have  been  taken  advantage  of  in  the  present  investiga¬ 
tion.  Where  the  integration  is  from  -<*>to  ♦  <*>,  all  odd  terms* 

drop  out.  Thus,  for  theM-2  case  .calculations  for  and 
result  in  the  same  3x3  matrix.  In  all  the  analyses  in  this 
report,  therefore,  the  eigenvalue  matrices  are  psuedosymmetric ; 
that  is,  they  are  symmetric  if  one  considers  that  the  3x3 
submatrices  are  the  eigenvalue  matrix  elements. 

This  set  of  t\  calculations  is  carried  out  by  the  computer 
program  in  terms  of  n  values  of  the  parameter  K,  as  shown  in 
Table  II;  the  value  of  K  corresponds  to  the  number  of  matrix 
elements  that  must  be  calculated. 


*  Here  "odd  terms"  refer  to  those  which  change  sign  in 
crossing  zero  on  the  axis  of  integration. 
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TABLE  I.  VALUES  OF  -  MU  -  %  ,  THE  PHASE 

DIFFERENCE  CONSTANT  FOR  n-BLADED  ROTOP. 
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TABLE  II.  NO.  OF  MATRIX  ELEMENTS  RESULTING  FROM  EACH  CALCULATION  K 

‘Lj, 

No.  of  matrix  elements  re¬ 
sulting  from  this  set  of 
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3.4  EIGENVALUE  ANALYSIS 


For  a  set  of  n  vortex  helices,  the  perturbation  equations  are  a 
set  of  nxn-  self-  and  mutual -inductance  effects.  Each  of 
these  nxn  components  in  turn  is  a  3x3  matrix  involving  per¬ 
turbations  in  the  r ,  <£  and*  directions.  The  resulting  per¬ 
turbation  equations,  therefore,  constitute  a  3nX  3n  eigen- 
matrix,  which  must  be  solved  to  conduct  the  stability  analysis. 
Each  3x3  matrix  has  complex  elements.  As  suggested  by  use  of 
the  term  "psuedo-symmetric"  in  the  previous  section,  although 
the  n  xn  array  of  matrix  elements  is  symmetric,  when  the  3x3 
matrix  elements  are  inserted,  the  3nx  3n  stability  matrix 
is  not  symmetric. 

The  eigenmatrix  equation  which  expresses  the  conditions  for 
solution  of  admissible  perturbation  growth  rates  can  be  written 
as 

A  x  =  Ax 

where  *  ■  3nxi  vector  corresponding  to  the  3n  perturbation 

displacement  components. 
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A  ■  nonsymmetric  3nx3n  matrix  with  complex  elements. 

A  ■  eigenvalues  of  the  above  matrix  equation;  a  set  of 
3n  complex  quantities,  related  to  the  dimensional 
growth  rates  and  wave  velocities  of  the  perturba-  a 
tions,  i.e.,  the  complex  quantity,  o<.  ,  by  >  -  . 

Solving  this  perturbation  equation  is  a  substantial  task  even 
with  the  help  of  a  digital  computer.  &  is  nonsymmetric,  has 
complex  elements, and  has  the  order  18x18  for  a  six-bladed 
rotor. 

Considerable  effort  was  expended  in  attempts  to  simplify  the 
above  matrix  equation.  Section  2.8  describes  one  such  simpli¬ 
fication  for  a  two-bladed  rotor.  This,  however,  cannot  be  gen¬ 
eralized  for  three-,  four-,  five-,  and  six-bladed  rotors,  since 
the  structure  of  the  eigenvalue  matrix,  £  ,  does  not  permit 
a  transformation  of  the  type 

*•!  ~ 

where  ^  ■  a  new  perturbation  vector  (  for  the 

vortex. 

X*  ■  a  perturbation  vector  (  Stv,y&4v  for  the 

vortex. 

a  set  of  scalars. 

The  subroutine  labled  ' ALLMAT '  has  been  used  to  determine 
eigenvalues  and  eigenvectors.  ’ ALLMAT '  is  a  FORTRAN  IV  subpro¬ 
gram  that  calculates  the  right  eigenvalues*  and/or  eigenvectors 
of  arbitrary  complex  matrices,  using  the  QR  algorithm  and  the 
Wielandt  inverse  power  method  for  vectors. 

Some  discussion  concerning  the  suitability  of  this  subroutine, 
the  technique  used  for  calculations,  and  its  numerical  stability 
seems  pertinent.  The  theoretical  background  for  the  QR  algo¬ 
rithm  may  be'  found  in  numerical  analysis  texts,  for  example, 
Ralston  [Reference  24  J  and  Acton  [Reference  25],  The  function 
of  the  QR  algorithm  is  to  decompose  an  arbitrary,  nonsymmetric 
complex  matrix  A  into  a  product  ay,  where  $  is  unitary  and 
u  is  upper  triangular.  Since  this  method  utilizes  unitary 


*  IfAis  a  nonsymmetric  complex  matrix  and  4*  denotes  its  Her- 
mitian  transpose,  right  eigenvectors  and  left  eigenvec¬ 
tors  jjj  are  deiined  by  relations  An  and  y?A  =  v,y!H 

respectively.  "  “ 
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transformations,  as  opposed  to  triangular  decomposition,  it 
tends  to  be  numerically  stable.  An  alternate  approach  to 
treating  the  complex  matrix  eigenvalue  problem  is  to  reduce  the 
problem  to  one  involving  real  matrices.  This  technique, 
however,  doubles  the  order  of  the  matrix,  and  for  large-order 
matrices  it  is  generally  preferable  to  treat  a  complex  matrix 
directly. 

The  complex  subroutine  'ALLMAT '  calculates  the  full  set  of 
eigenvalues  (3rt  in  number),  after  the  eigenmatrix  elements 
have  been  evaluated  for  the  specified  values  of  the  parameters 
fc  (pitch),  a>  (wave  number  of  the  perturbations)  and  fc  (vortex 
core  size  specified  as  a  fraction  of  helix  diameter  or  radius) 
for  an  n-bladed  rotor.  These  eigenvalues,  in  general,  are 
complex.  The  real  part  of  the  eigenvalues  determines  the  decay 
or  the  growth  rates  of  the  perturbations,  thereby  showing 
stability  or  instability,  respectively.  Positive  real  eigen¬ 
values  are,  of  course,  unstable  while  negative  real  eigenvalues 
are  stable.  Eigenvalues  with  zero  real  parts  are  "neutrally 
stable".  It  is  sufficient,  therefore,  to  examine  the  eigen¬ 
value  with  maximum  real  part  to  determine  stability.  The 
imaginary  parts  of  the  eigenvalues  define  the  direction  of  prop¬ 
agation  of  the  disturbances.  An  eigenvalue  with  a  negative 
imaginary  part  corresponds  to  disturbances  propagating  in  the 
positive  ©  direction. 

The  stability  analysis  results  obtained  in  this  study  indicate 
that  there  are  in  all  cases  disturbances  travelling  in  several 
directions  with  positive,  negative  and  zero  growth  rates.  In  a 
linear  stability  analysis  such  as  this,  the  unstable  modes  with 
the  largest  divergence  rates  very  quickly  dominate  the  distur¬ 
bance  shapes  and  determine  the  shape  of  the  resulting  distorted 
helices. 


CHAPTER  4 


RESULTS 


4.1  GENERAL 


This  chapter  presents  three  distinct  kinds  of  material.  The 
first  consists  of  a  discussion  of  the  selection  of  cases  and  a 
physical  description  of  the  modal  deflections.  These  are 
contained  in  Sections  4.2  and  4.3.  The  second 
consists  of  the  results  for  special  cases  which  were  selected 
to  help  verify  the  theory  and  increase  confidence  in  the 
numerical  methods  used  in  the  present  investigation.  This 
material  is  contained  in  Sections  4.4  through  4.6.  The  third 
kind  of  material  is  comprised  of  results  for  two-,  three-, 
four-,  five-,  and  six-bladed  rotors  for  various  specific  values 
of  parameters  of  interest,  such  as  the  pitch  of  the  helix, 
the  wave  number  of  perturbations,  and  core  size.  These  results 
are  presented  in  Sections  4.7  and  4.8  . 

There  is  very  little  experimental  evidence  against  which  to 
check  the  predictions  carried  out  in  the  present  investigation 
of  wake  stability.  However,  some  qualitative  experimental 
results  have  been  reported  by  Landgrebe  [Reference  11]  using  a 
rotor  model  hover  test  facility.  The  object  of  Landgrebe' s 
experiment  was  to  measure  the  effect  of  parameters  such  as 
number  of  blades,  blade  linear  twist,  blade  aspect  ratio, 
rotor  tip  speed  and  collective  pitch  on  rotor  hover  performance 
and  the  associated  wake  geometry  characteristics.  Flow  visual¬ 
ization  data  were  obtained  by  injecting  smoke,  and  wake 
photographs  and  high- frame- speed  movies  were  taken.  Examination 
of  the  photographic  results  leads  to  the  qualitative  conclusion 
that  the  far-wake  region  of  a  hovering  rotor  is  unstable  or, 
at  best,  neutrally  stable.  These  conclusions  do  not,  at  least, 
contradict  the  results  of  the  present  investigation. 


4.2  SELECTION  OF  CASES  TO  BE  RUN 


In  this  section  representative  values  of  the  pitch  of  the 
helix  have  been  calculated  under  hovering  conditions.  Glauert, 
for  example  [Reference  1],  and  Payne  [Reference  26]  give  the 
following  equation  for  the  induced  velocity  at  the  rotor  in 
vertical  flight. 

Wc+'\)'>L=  ^ 


where  vc  *  climb  velocity. 
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■  increase  in  velocity  at  the  actuator  disc  (  induced 
velocity) . 

2%.  ■  disc  loading. 

T 

P  ■  mass  density  of  air  ■  0.0023  slug/ft  (sea  level, 
std.  day). 

In  hover,  Vc  «  0;  therefore i  vt  -  far.  The  downward  axial 

\/2fi 

velocity  in  a  fully  developed  slipstream  V*  is  given  by 


Vi  -  21 >t 


2 

Typical  helicopter  disc  loadings  vary  from  2  to  10  lb/ft  , 
and  typical  tip  speeds  range  from  500  to  700  ft/sec. 

Thus,  since  average  axial  downwash  velocity  is  given  by 


typical  values  vary  from  40  ft/sec  to  95  ft/sec, 

A  period  for  one  rotor  revolution  is  T  ,  and  axial  dis 
placement  for  one  turn  of  the  helix  is  f  -hre  ■  krun)  , 
Therefore  , 

~  V,  •  kr  irr 

°T  k  m  It  m  axial  downwash  speed 

Jlr  tip  speed  of  roror 


Thus  values  of  pitch  of  the  wake  helix  corresponding  to 
typical  disc  loadings  and  tip  speeds  will  vary  from 


70$  B  °*os  t0  3J0*  0-20 
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Although,  as  shown  above,  the  pitch  of  the  vortex  helices 
corresponding  to  the  normal  extremes  of  hovering  rotor  opera- 
ting  conditions  lies  approximately  between  0.05  and  0.20,  the 
majority  of  results  obtained  here  are  for  pitch  equal  to  0.1 
radian.  The  ratio  of  vortex  core  diameter  to  the  helix  diam¬ 
eter,  labeled  in  the  tables  given  in  this  report  as  "  core 
size", is,  however,  considerably  more  difficult  to  relate  to 
reality.  A  paucity  of  experimental  results  or  other  evidence 
exists  regarding  actual  vortex  core  size  of  the  helices  in  the 
wake  of  a  hovering  rotor.  The  ratio  of  core  radius  to  helix 
radius  has  been  chosen  here  to  be  0.1,  following  Widnall,  and 
some  cases  have  also  been  for  this  nondimensional  core  size 
equal  to  0.33  to  investigate  the  effect  of  core  size  on  vortex 
stability.  Perturbation  wave  numbers  have  been  examined  from 
0.25  (l<mg  wavelengths)  to  8.0  (short  wavelengths),  including 
the  wave  number  of  0.0.  Perturbation  wavelengths  which  are 
very  short  compared  to  the  helix  diameter  are  likely  to  occur 
if  the  disturbances  are  concentrated  and  would  seem  to  involve 
changes  of  vortex  core  cross  section  in  the  plane  perpendicular 
to  the  helical  filament  centerline.  Thomson  [Reference  27], 
however,  has  shown  that  so  long  as  the  radius  of  the  cross 
section  of  the  core  in  the  plane  perpendicular  to  the  filament 
centerline  is  small  compared  with  the  radius  of  curvature 
of  that  centerline,  the  cross  section  will  remain  approxi¬ 
mately  circular  under  the  self- influence  of  the  whole  vortex. 
This  same  observation  has  been  made  by  Widnall  in  her  study  of 
short  wave  perturbation  instability  modes  of  planar  circular 
vortex  rings  [Reference  28];  however,  a  cutoff  at  rather  large 
wave  numbers  was  used  to  limit  the  extent  of  the  stability 
curves  in  Reference  16.  The  core  diameter  ratio  of  0.10 
satisfies  the  assumption  that  the  core  size  is  small  compared 
to  helix  radius;  0.33  does  not.  For  the  highest  wave  numbers, 
i.e.,  8.0,  even  a  core  which  is  0.1  helix  diameter  is  hardly 
small  compared  to  a  wavelength  of  0.393  helix  diameter.  Thus, 
in  Figures  14  through  19  a  wave  number  limit  of  about  8.0  is 
arbitrarily  shown. 


4.3  PHYSICAL  DESCRIPTION  OF  MODAL  DEFLECTIONS 


A  clear  picture  of  modal  perturbations  will  provide  increased 
insight  into  instability  mechanisms.  The  type  of  vector 
perturbations  introduced  in  Section  2.6  is  represented  by 


«t+  lie  6m 


This  deflection,  of  course,  has  components  in  r  ,  <Jj  and*  di¬ 
rections.  To  visualize  the  deflections  associated  with  certain 
wave  numbers,  consider  a  deformation  of  the  type  &r  *  Srcosiuei 
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where  6v  *  perturbation  in  r  direction. 

A 

fcr  ■  amplitude  of  the  perturbation, 
to  ■  the  wave  number. 

The  real  parts  of  distortion  modes  in  axial  and  circumferential 
directions  can  similarly  be  written  as 

s4>  =  ICosco©] 


and 


Si  =  Si  Cco-S  to©]  e. 


Kfc 


Once  the  eigenfunctions  and  Si  are  known,  the 
distorted  shapes  of  the  helices  in  these  directions  can  be 
calculated  for  the  corresponding  wave  number, to  . 


Dilatational  Mode,  03  -  0 
For  to  ■  0  , 

V  =  6?  e“‘ 

A 

scf>  =  s<£  e 
62  =  62:  e 

Since  this  perturbation  is  independent  of  6,  it  is  uniform 
throughout  space.  Whereas  the  i  and  4>  coordinates  simply 
change  the  position  of  given  helix  with  respect  to  coordinate 
axes,  the  radial  perturbations  expand  or  contract  the  diameter 
of  the  vortex  helix  uniformly.  The  unperturbed  helix  is 
shown  in  Figure  1,  and  a  dilatational  mode  disturbance  is 
shown  in  Figure  7.  Note  that  to  simplify  sketches  of  distor¬ 
ted  wake  shapes,  in  all  cases,  perturbations  will  be  shown  in 
only  one  direction  in  any  one  sketch.  This  can  be  done  without 
sacrificing  generality,  so  long  as  the  reader  keeps  in  mind 
that  reference  to  more  than  one  figure  will  generally  be  re¬ 
quired  to  visualize  a  given  modal  deflection.  Since  this  dila¬ 
tational  mode  distorts  the  helix  from  its  initial  configuration 
only  by  uniform  diameter  changes,  it  is  a  disturbance  which, 
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in  a  nondimensional  sense,  changes  the  helix  pitch  angle 
without  changing  the  basic  helica]  shape. 


Deformation  Mode  for  ■  1 

For  co  m  if  the  deformation  shape  can  he  pictured  by  noting  that 
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These  deformations  are  sketched  in  Figures  8  and  9. 

As  shown  in  the  top  and  side  views,  there  is  some  deformation 
of  the  initial  circular  shape  of  the  helix  as  a  result  of 
perturbations  in  y  direction,  but  in  the  main,  it  represents 
a  relatively  undistorted  translation  in  a  direction  normal  to 
the  centerline  axis  of  the  helix.  In  the  e-  direction,  however, 
these  deformations  amount  to  a  tilting  of  the  coils  relative 
to  the  axis  of  the  helix.  This  point  has  been  elaborated 
in  Section  4.4.  Deformations  in  the  <P  or  circumferential  sense 
mainly  stretch  or  compress  the  helical  vortex  line,  and  since 
in  this  analysis  the  cross  section  of  the  vortex  core  (i.e., 
normal  to  the  helical  axis)  has  been  assumed  to  be  constant, 
deformations  in  the  circumferential  sense  can  largely  be 
ignored. 


Distortion  Mode  Values  of  oo  Other  Than  0  or  1 
Consider,  for  example,^"  2 
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The  plan  and  side  views  of  the  helix  distorted  by  radial  and 
axial  deformations  in  this  mode  are  shown, respectively ,  in 
Figures  10  and  11.  In  Figure  10,  the  helix  is  shown  to  be 
distorted  from  circular  to  elliptic.  In  Figure  11,  the 
distortions  change  the  projections  of  the  helical  turns  from 
straight  lines  to  sickle-shaped  curves,  which  are  all  either 
concave  upward  or  concave  downward,  depending  upon  the  sign  of 
the  deformation  in  the  -a-  direction. 


Other  Possible  Distortion  Modes 


Distortion  forms  corresponding  to  higher  integer  and  noninte¬ 
ger  wave  numbers  can  be  similarly  visualized.  For  example,  the 
distorted  shapes  of  two  successive  coils  are  quite  different 
when  to  ■  1/2,  but  such  distortions  are  repeated  after  traver¬ 
sing  every  two  coils  of  the  helix.  Similarly,  for  uo  »  1.5, 
three  coils  of  the  helix  must  be  examined  before  one  would  ex¬ 
pect  to  find  a  complete  cycle.  On  the  other  hand,  for  co  ■  3, 
there  will  be  three  maximum  and  three  minimum  distortion  points 
in  one  turn  of  helix.  The  co  ■  1/2  mode  for  r  and  i  pertur¬ 
bations  is  shown  in  Figures  12  and  13a. 

Extension  of  Modal  Perturbation  Modes  to  Multiple  Interdigita- 
ted  Helices 


Each  vortex  helix  in  the  array  trailed  by  multibladed  rotors 
can  undergo  the  modal  deformations  discussed  in  the  previous 
sections.  In  addition,  however,  the  temporal  phase  relation¬ 
ship  among  the  deformations  in  different  helices  provides 
three  more  possibilities  of  combinations  of  motions,  one  for 
each  coordinate  direction,  for  each  additional  helix.  The 
actual  phase  and  magnitude  are  determined  by  the  eigen-solution. 
In  the  simplest  possible  example,  two  interdigitated  helices 
may  be  undergoing  a  dilatational  (i.e.,  0)  deformation,  but 

with  opposite  sign,  so  that  when  one  helix  is  expanding  the 
other  is  contracting, and  vice  versa.  Similarly,  in  a  no- 
radial-deformation  (i.e.,co»  1)  mode,  the  modal  deformations 
for  each  helix  of  a  two-bladed  rotor  may  have  the  same  sign, 
in  which  case  the  two  helices  would  translate  in  the  same 
direction,  and  the  change  in  pitch  angle  would  also  be  occur¬ 
ring  in  the  same  direction.  In  this  case  there  would  be  no 
change  in  the  relative  positions  of  the  two  helices,  one  to  the 
other.  By  contrast,  if ,  in  an co  =  l  mode,  the  modal  deforma¬ 
tions  of  the  two  helices  had  opposite  signs  and  equal  magni¬ 
tude,  then  the  translation  of  each  helix  would  be  in  opposite 
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directions  and  the  effective  pitch  change  of  these  helices 
would  also  be  such  as  to  reduce  the  axial  distance  between 
corresponding  segments  of  the  helices. 


4.4  LEVY  AND  FORSDYKE  CASE 


The  early  Levy  and  Forsdyke  analysis  [Reference  15]  cited  ear¬ 
lier  in  this  report  concluded  that  a  doubly  infinite,  constant- 
diameter  helical  vortex  is  unstable  for  the  wave  number  ■  1  mode 
when  the  helix  pitch  angle  is  less  than  0.3  radian  and  is 
stable  for  pitch  angles  greater  than  0.3  radian.  Stability 
analyses  for  the  same  single  helix  case  have  also  been  carried 
out  in  the  present  study  and  have  shown  such  a  vortex  to  be 
stable  in  the  wave  number  *  1  perturbation  mode  for  all  helix 
pitch  angles.  Some  insight  is  gained  by  reference  to  the  dis¬ 
tortion  shapes  in  this  mode. 


The  *1  perturbation  comes  close  to  translating  the  whole 
helix  to  right,  left.upior  down  without  much  radial  distortion. 
Such  motion  clearly  is  approximated  by  a  change  in  the 
arbitrary  choice  of  a  coordinate  axis  system,  and  it  cannot 
reasonably  be  expected  to  influence  the  tendency  of  the 
helix  to  distort.  Such  motion  is  depicted  in  Figure  8. 

However,  deformations  in  the  direction  presented  in  Figure 
9  for  this  mode  are  shown  to  be  equivalent  to  tilting  the 
coils  relative  to  the  axis  of  the  helix.  In  this  distortion, 
segments  at  any  point  on  the  helix  always  remain  at  the  same 
distance  from  those  segments  in  preceding  or  subsequent  turns 
directly  above  and  below  them.  With  respect  to  their 
separation  from  vortex  segments  at  different  azimuthal  stations 
on  the  helix,  any  given  segment  gets  a  bit  closer  to  those 
on  one  "turn"  and  that  same  amount  farther  away  from  the 
corresponding  segment  on  the  next  "turn".  Thus,  it  is 
intuitively  satisfying  to  find  that  such  deformation  modes  are 
predicted  to  be  stable. 

Secondly,  beyond  those  discrepancies  in  the  Levy  and  Forsdyke 
results  associated  with  the  use  of  a  planimeter,  referred  to 
earlier  in  this  report,  their  evaluation  of  singular  integrals 
expanded  Cosw  and  terms  in  power  series.  Terms  up  to 

*  were  retained  in  the  numerators  and  terms  up  to  x  were  re¬ 
tained  in  the  denominators.  Since  the  upper  and  lower  limits 
of  integration  are  0  to  */(>  and  the  integrands  are  near  singu¬ 
lar  in  their  behavior,  considerable  error  is  introduced  by  this 
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approximation.  Table  III  lists  cases  run  under  the  present 
investigation  which  may  be  compared  with  the  results  of  Levy 
and  Forsdyke  [Reference  15].  In  addition,  selected  integrals 
were  calculated;  most  of  the  integrated  values  checked  out, 
but  errors  as  high  as  an  order  of  magnitude  were  noted  in 
others. 


4.5  WIDNALL  CASE  OF  SINGLE  HELIX 


The  results  recently  published  by  Widnall  [Reference  16]  repre¬ 
sent  the  stability  of  a  single  vortex  filament  of  finite  core 
for  small  sinusoidal  displacements  of  its  (helical)  center- 
line.  The  singularity  resulting  from  the  self- induced  motion 
of  the  element  has  been  removed  in  Reference  16  by  cutting  off 
the  integration  in  the  region  of  the  singularity  and  using  a 
method  of  matched  asymtotic  expansions  in  this  region.  The 
solution  in  this  region  has  been  matched  with  the  solution  for 
a  plane,  circular  vortex  ring  of  the  same  radius  of  curvature. 
This  mathematically  sophisticated  approach  succesfully  dupli¬ 
cated  the  results  of  Betchov  [Reference  29]  for  the  case  of  a 
vortex  filament  whose  cross-sectional  radius  is  very  small 
compared  to  its  radius  of  curvature.  Numerical  results  for 
various  pitch  angles,  core  sizes  and  perturbation  wave  numbers 
have  been  reported  by  Widnall. 

The  special  case  of  a  single  helix  in  the  present  study  deals 
with  the  same  physical  situation  as  the  work  reported  in  Ref¬ 
erence  16  but  differs  significantly  in  the  treatment  of 
singularities,  as  has  been  mentioned  in  chapter  2. 

The  results  of  the  single  helix  cases  run  under  the  present 
study  are  plotted  in  Figure  14  and  tabulated  in  Table  IV 
along  with  the  corresponding  results  obtained  by  Widnall.  The 
table  facilitates  a  direct  comparison  of  the  numerical  results 
of  this . analysis  with  those  of  Widnall.  One  observes  that  the 
results  are  qualitatively  similar,  but  that  the  numerical 
values  of  maximum  divergence  rates  differ.  In  a  typical  case, 
for  example,  pitch  *  0,1,  wave  number  =0.5  and  core  size  ■ 
0.33,  Widnall  predicts  a  divergence  rate  of  24.8  The  corres¬ 
ponding  value  in  the  present  analysis  is  7.90.  Note,  however, 
that  in  the  present  investigation  divergence  rate  is  defined  as 
,  while  in  the  Widnall  study  divergence  rate  is  taken  to 
be  Thus,  for  a  comparison,  the  present  value  of  7,90 

must  be  compared  with  (1/2)24.8  =  12.4.  Thus,  to  faci¬ 
litate  comparisons  on  the  same  nondimensionalized  basis,  the 
divergence  rates  taken  from  the  Widnall  graphs  have  been  divided 
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TABLE  III.  LEVY  AND  FORSDYKE  CASE,  COMPARISON 

WITH  THE  PRESENT  STUDY 

■  i 

Stable  or 

Unstable 

Levy/ Forsdyke 

Present 

Wave  Number 

Core  Size 

Pitch 

Study 

Study 

1.00 

0.10 

0.10 

Unstable 

Stable 

1.00 

0.10 

0.20 

Unstable 

Stable 

1.00 

0.10 

0. 25 

Unstable 

Stable 

1.00 

0.10 

0.30 

Stable  or 
Unstable 

Stable 

1.00 

0.10 

0.40 

Stable 

Stable 

1.00 

0.10 

0.60 

Stable 

Stable 

1.00 

0.10 

0.80 

Stable 

Stable 

1.00 

0.10 

1.00 

Stable 

Stable 

'  . -  -  —  - .  .  _  —  i 
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TABLE  IV.  CASE  OF  SINGLE -BLADE ,  DIVERGENCE  RATES 
OBTAINED  IN  THE  PRESENT  STUDY 
COMPARED  WITH  THE  VALUES  OBTAINED 
BY  WIDNALL  (Ref.  16*) 


Wave 

Number 

Core  Siz< 

8  ■  0 

.1 

Max.  Diver¬ 
gence  Rates 
in  Present 
Study 

Max.  Diver¬ 
gence  Rates 
Reported  by 
Widnall 

0.00 

-0.00 

0.0 

(0.0) 

0.25 

8.75 

19.0 

(9.5) 

o.so 

11.84 

25.2 

(12.6) 

0.75 

9.01 

18.8 

(9.4) 

1.00 

0.00 

0.0 

(0.0) 

1.25 

8.33 

18.0 

(9.0) 

1.50 

11.43 

24.0 

(12.0) 

1.75 

8.80 

16.6 

(8.3) 

2.00 

0.00 

0.0 

(0.0) 

2.50 

10.61 

21.4 

(10.7) 

3.00 

0.26 

0.0 

(0.0) 

3.50 

9.47 

0.0 

(0.0) 

4.00 

0.99 

0.0 

(0.0) 

4.50 

8.23 

0.0 

(0.0) 

5.00 

1.99 

0.0 

(0.0) 

5.50 

7.18 

0.0 

(0.0) 

6.00 

3.21 

0.0 

(0.0) 

7.00 

4.51 

0.0 

(0.0) 

8.00 

5.79 

0.0 

(0.0) 

0.0 

17.6 
24.8 

17.4 
0.0 

16.4 
23.2 

16.4 
0.0 

22.6 
7.6 

20.6 

12.0 

19.6 
14.0 

17.4 

11.6 
0.0 
0.0 


(0.0) 

(8.8) 

(12.4) 

(8.7) 
(0.0) 
(8.2) 
(11.6) 
(8.2) 
(0.0) 

(11.3) 

(3.8) 

(10.3) 

(6.0) 

(9.8) 
(7.0) 

(8.7) 

(5.8) 
(0.0) 
(0.0) 


*  The  approximate  divergence  rates  used  in  the  above  com¬ 
parison  were  taken  from  the  graphs  in  Ref.  16, 
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by  a  factor  of  2  and  written  within  the  brackets  in  the  same 
columns  in  Table  IV.  As  pointed  out  in  the  previous  paragraph, 
the  results  are  significantly  different  for  core  size  0.33;  the 
divergence  rates  compare  quite  closely,  however,  for  core 
size  0.1.  It  has  been  noted  that  estimates  of  the  core  size  of 
vortex  helices  in  the  wake  of  a  hovering  rotor  are  not  known 
with  any  confidence,  but  a  core  diameter  ratio  of  0.1  would 
seem  closer  to  reality  than  0.33.  Note  that  Landgrebe  [Refer¬ 
ence  11]  assumed  the  core  radius  of  0.005,  while  Crimi  [Refer¬ 
ence  12]  has  used  a  value  of  0.05  to  0.08 

Some  estimates  of  core  size  are,  however,  available  for  the 
vortices  trailed  by  fixed-wing  aircraft.  Spreiter  and  Sacks 
[Reference  30]  estimated  the  core  diameter  of  vortices  trailed 
behind  elliptically  loaded  aircraft  wings  as  0.197  times  the 
separation  between  the  trailing  vortices.  The  effect  of  curva¬ 
ture  probably  modifies  the  core  size,  but  this  estimate  of  core 
size  falls  somewhat  closer  to  the  0.1  value  used  for  helical 
vortices  in  the  present  study  than  it  does  to  0.33.  For  the 
smaller,  more  realistic  core  sizes,  the  present  study 
with  its  different  approach  toward  treatment  of  singularities 
in  the  Biot-Savart  integration,  appears  to  predict  divergence 
rates  quite  comparable  to  the  method  used  in  the  Widnall  study. 


4.6  S.  C.  CROW’S  CASE  (TRAILING  VORTICES  FROM  FIXED-WING 

- Aiirnr  wTNcsy —  - 


S.  C.  Crow  [Reference  17]  presents  an  analysis  of  the  stability 
of  a  pair  of  straight,  parallel  vortices  as  trailed  from  a  fixed- 
wing  aircraft.  This  analysis  also  gives  rise  to  an  eigenvalue 
problem  for  the  growth  rate  of  sinusoidal  perturbations. 

If  the  pitch  is  allowed  to  approach  infinity  in  the  present  stu¬ 
dy  for  the  case  of  a  two-bladed  rotor,  the  wake  should  correspond 
to  two  straight,  parallel  infinite  vortices.  It  is  shown  in 
Appendix  I  that  the  perturbation  equations  used  here,  taken  to 
this  limit,  yield  the  same  expressions  as  S.  C.  Crow's  case,  ex¬ 
cept  for  tne  differences  associated  with  the  cutoff  technique 
to  remove  singularities.  Incidentally,  P.  C.  Parks  [Reference 
20]  has  duplicated  S.  C.  Crow's  results  by  modifying  the  self- 
induction  integrals  in  lieu  of  "cutting  off".  Parks'  integra¬ 
tion  to  eliminate  the  effects  of  singular  integrals  yields  the 
same  expressions  obtained  in  the  present  study,  after  the  proper 
limits  are  taken,  and  has  confirmed  the  essential  features  of 
Crow's  theory  with  small  numerical  changes.  Since  the  singular¬ 
ities  in  the  self-induction  integral  have  been  removed  in  the 
present  study  by  an  approach  similar  to  that  used  by  P.  C. 

Parks,  it  is  reassuring  that  limiting  expressions  which  are 
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approaching  the  case  of  trailing  fixed-wing  vortices  do,  in 
fact,  duplicate  those  used  by  Parks. 

It  is  a  more  difficult  matter  to  attempt  to  duplicate  S.  C. 
Crow's  quantitative  results  using  the  computer  program  devel¬ 
oped  in  the  present  study,  since  the  integration  techniques 
used  here  are  tailored  for  low  values  of  pitch.  In  limiting 
cases  where  pitch  is  taken  to  be  very  large,  the  sine  and 
cosine  terms  in  the  perturbation  equations  oscillate  very 
rapidly, and  the  integrated  values  using  the  current  numerical 
techniques  are  inappropriate  for  obtaining  accurate  results. 

Two  cases,  however,  have  been  run, and  the  results  are  tabulated 
in  Table  V. 


4.7  DISCUSSION  OF  RESULTS 


Single  Helix 

The  single-blade  case  was  run  primarily  to  check  the  present 
work  against  the  results  of  Levy  and  Forsdyke  [Reference  15] 
and  Widnall  [Reference  16],  The  comparisons  have 
been  discussed  in  the  preceding  sections.  Further  considera¬ 
tion  of  the  single  helix  case  will  be  limited  here  to  a 
discussion  of  the  perturbation  deflection  shapes  and  their 
influence  on  helix  stability.  As  shown  in  Figure  14,  the 
single  helix  has  no  instabilities  for  wave  numbers  «*>*  0  and 
w  ■  1.  Figures  7,  8,  9  and  13d  show  that  radial  and  axial 
deflections  for  these  two  deformation  modes  have 
two  characteristics  in  common:  first,  the  circular  cross 
section  seen  in  the  plan  views  is  largely  preserved,  and 
second,  axial  separation  between  vortex  segments  at  correspon¬ 
ding  circumferential  positions  is  preserved  in  the  deformations. 
In  terms  of  the  Biot-Savart  law  of  induction  it  is  clear  that 
changes  of  the  helix  diameter,  so  long  as  they  are  uniform 
around  the  circumference, will  cause  uniform  changes  of  induced 
effects.  Further,  if  such  deformations  are  repeated  without 
change  on  every  helix  coil,  then  induced  effects  associated 
with  changes  in  helix  diameter  will  not  distort  the  helix,  but 
will  simply  cause  the  entire  helix  to  translate  axially.  Since 
the  choice  of  coordinate  axes  and  their  velocities  in  inertial 
space  is  completely  arbitrary,  it  is  not  surprising  that 
neither  translation  of  the  entire  helix,  such  as  exists  in  the 
plan  view  for  co  ■  l ,  nor  uniform  change  of  the  helix  radius,  as 
shown  in  the  plan  view  of  Figure  7  for  to  =  0,  has  any  in¬ 
fluence  on  the  stability.  Changes  in  the  relative  axial  posi¬ 
tion  of  helix  segments,  however,  will  influence  induced  effects 
in  the  radial  direction.  For  these  two  modes,  however,  there 
is  no  such  change  in  relative  position,  at  least  among  the 
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TABLE  V.  RESULTS  OF  THIS  ANALYSIS  FOR  S.  C.  CROW'S  CASE 


L _ 


Point  1 

Point  2 

Eigenmatrix 

Eigenmatrix 

0 

3.92263 

0 

-0.91889 

0 

D. 85992 

0 

-3.5407 

0 

0 

0 

-1.7587 

O 

0 

►0.91889 

0 

0 

3.85935 

0 

-0.85992 

♦1.6991 

0 

0 

0 

+3.4780 

0 

A 

*i*[(-0.qii6  3  +  0.'n«9J(J.5*o7-l.6Qqi)J  1 

*  Stable. 

■[_<- 0 . q »  - 0 .  <7 fc-4-0?  - 1  6  9 9  1  )  3  ''x 

=  C-\»e3'/x  Stable. 

y 

*s*U-o.9s<m.t-o,95qssxs'4>8o-iM«?)lJ 
so  Neutrally  Stable. 

*s=  tc-o.96‘?<u-o.«sq3sXa.478o+i-«*?)]  2 

i j 

=  C-v«J  x  Stable. 

Point  1: 

Parameters  in  S.  C.  Crow' 

s  paper: 

§  -2.  %-°-1 

Point  2: 

Parameters:  Q  *0.5,  5/^, 

“0.1 

corresponding  parameters  in  present 

study 

Point  1: 

Pitch  «  1000,  wave  number 

=  1000, 

&  *=  0.311 

Point  2: 

Pitch  *  1000,  wave  number 

=  250  , 

e  *  0.311 

Remarks:  Numerical  comparison  of  divergence  rates  with  Crow's 
results  for  these  points  has  not  been  shown  here 
since  the  current  program  is  not  designed  for  such 
large  values  of  pitch  and  wave  numbers;  as  a  result, 
the  calculated  values  are  not  expected  to  be 
accurate.  The  purpose  of  the  above  run  is  to  demon¬ 
strate  the  limiting  behavior  in  the  present  study, 
which  has  otherwise  been  shown  theoretically. 
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vortex  segments  closest  to  one  another,  that  is,  in  the  same 
circumferential  position  from  one  coil  to  the  next.  For 
w  ■  1/2. on  the  other  hand,  as  shown  in  the  side  view  of 
Figure  12,  the  relative  axial  displacements  of  vortex  segments 
at  the  same  circumferential  position  do,  in  fact,  approach 
each  other  or  separate  from  one  another.  Induced  effects  of 
some  importance  do,  therefore,  exist,  and  any  initial  tendency 
for  adjacent  segments  to  come  together  will  be  amplified  still 
more  by  the  increased  proximity.  One  would  certainly  expect 
such  modes  to  be  unstable,  and  they  are,  as  shown  in  Figure  14. 

For  the00*  2  mode,  the  plan  view  changes  as  shown  in  Figure  10 
from  circle  to  ellipse.  Such  changes  induce  axial  velocities 
which  are  not  uniform  around  the  circumference.  Furthermore, 
they  will  be  such  as  to  cause  roughly  the  same  displacement 
changes  as  seen  in  the  side  view  of  modal  perturbations  shown 
in  Figure  10.  These  are  not  very  strong  effects,  however, 
since  for  every  pair  of  opposite  segments  which  come  together 
in  the  sense  seen  in  the  plan  view  (and  thereby  increase  their 
mutually  induced  effects) , there  are  two  other  segments  which 
are  farther  apart, not  just  from  each  other,  but  from  the  first 
pair.  Since  the  induced  effects  of  the  latter  pair  on  the 
first  pair  are  reduced,  there  are  tendencies  toward  compensa¬ 
tion.  The  separation  distance  is  in  the  denominator  of  the 
Biot-Savart  law,  however,  so  that  such  effects  are  not  linear. 
This  possibly  explains  why  the  divergence  rate  in  the  co  *  2 
mode  is  less  stable  for  the  single  helix  whose  vortex  core 
diameter  ratio  is  .33  than  it  is  for  the  case  where  this 
parameter  is  .1,  as  shown  in  Figure  14. 

Finally,  Figure  11  shows  that  the  modal  perturbation  for co  ■  2 
preserves  the  axial  separation  between  points  at  corresponding 
circumferential  locations  on  the  single  helix. 


Multiple -Bladed  Rotors 

It  is  important  to  note  that  the  solution  method  for  multiple- 
bladed  rotors  is  identical  to  that  for  a  single  helix;  namely, 
a  wave  number,  co  ,  is  selected,  and  the  eigenmatrix  is  solved 
numerically  for  the  values  of  A  which  determine  the  divergence 
rates  corresponding  to  the  perturbation  modes  associated  with 
that  value  of  c o.  There  is,  of  course,  an  eigen  deformation 
mode  associated  with  each  eigenvalue,  A*  Components  of  tlje 
eigen  mode  for  the  single  helix  are  the  complex  values  of  6r  , 
S4>  and  .  Since  the  multi  bladed  rotor  introduces  addition 

al  degrees  of  freedom,  for  the  case  of  a  rotor  with  n  blades 
there  are  n  sets  of  eigenvectors,  srx*  ,  and  ,  where 

i  *  J,  2,  .  .  .  n,  for  every  eigenvalue,  A  .  The  phase 
relationship  between  motions  of  the  blades  i  =  1  and  i  *  2  is 
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embodied  in  the  complex  values,  say.  of  and  ;  this 
phase  relationship  will  determine  whether  the  helix  trailed 
from  the  number  1  blade  will  move  toward  or  away  from  the 
helix  trailed  from  the  number  2  blade  in  the  deformation  eigen* 
mode  corresponding  to  A  .  This  phase  relationship,  between 
motions  in  the  same  coordinate  direction  at  corresponding 
points  on  two  different  helices,  cannot  be  known  prior  to  solv¬ 
ing  the  eigenvalue  matrix.  This  contrasts  with  the  case  of  a 
single  helix,  for  which  phase  relations  among  segments  at 
corresponding  positions  are  known  as  soon  as  a  value  of  go  is 
chosen. 

Numerical  results  for  multiple -bladed  rotors  are  tabulated  in 
Tables  VI  through  X  and  are  plotted  in  Figures  15  through 
19.  Most  of  these  results  are  obtained  with  pitch  *  0.1  and 
core  size  ■  0.1.  A  few  cases  were  run  for  pitch  ■  0.15  and 
core  size  »  0.33;  these  are  also  shown  in  the  figures,  although 
the  curves  vs  wave  number  are  largely  extrapolated,  using  the 
curves  for  pitch  -  0.1  and  core  size  ■  0.1  as  guides. 

For  the  case  of  a  two-bladed  rotor,  the  divergence  rate,  X  ,  is 
plotted  vs  wave  number,  m  ,  in  Figure  15.  It  is  noted  that 
the  maxima  of  the  continuum  of  the  divergence  rate  points  de¬ 
crease  rather  smoothly  as  wave  number  increases.  Similarly  the 
minima  of  the  curves  increase  slightly  and  then  decrease  rather 
smoothly.  For  larger  integer  wave  numbers,  i.e.,  short  wave¬ 
length  perturbations,  the  difference  between  minimum  and 
maximum  divergence  rates  decreases  as  wave  number  increases. 
Figures  20  and  21  suggest  that  maximum  divergence  rates  occur 
at  integer  wave  numbers,  because  at  these  wave  numbers  the 
segments  adjacent  to  each  other  on  the  two  helices  approach 
each  other  to  the  maximum  extent.  Minimum  relative  motion  of 
this  sort  seems  to  occur  with  half -integer  wave  numbers. 
Imagining  a  helix  from  a  second  blade  in  Figure  13a  with  the 
most  favorable  phase  angle  possible  shows  how  little  of  the 
vortices  approach  each  other  in  such  motion. 

Figures  16  through  19.  for  three-,  four-,  five-,  and  six-bladed 
rotors,  respectively,  exhibit  similar  divergence  rate  -wave 
number  relationships;  i.e.,  the  maximum  and  minimum  divergence 
rates  occur  at  wave  numbers  equal  to  1/2  times  even-integer 
multiples  and  1/4  times  odd -integer  multiples,  respectively, 
times  the  number  of  blades.  For  example,  a  five-bladed  rotor 
exhibits  maximum  divergence  rates  at  the  wave  numbers  2.5,  5.0 
and  7.5,  etc.,  which  are  1/2,  1,  and  3/2  times  5.  Similarly, 
minimum  divergence  rates  for  a  five-bladed  rotor  occur  at 
wave  numbers  1.25,  3.75  and  6.25,  etc.,  which  is  1/4,  3/4 
and  5/4  times  the  number  of  blades  in  the  rotor .namely,  5. 

For  three-,  four-,  five- and  six-bladed  rotors,  the  excursion  from 
the  mean  divergence  rate  also  decreases  with  increasing  wave 
numbers. 
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TABLE  VI.  MAXIMUM  DIVERGENCE  RATES  FOR  TWO-BLADED  ROTOR 


Pitch  of 
Helix  (rad) 


Wave 

Number 


Max  Divergence  Rate 


46.65 

37.05 

46.26 

49.80 

46.57 

37.16 

45.42 

49.06 
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TABLE  VII.  MAXIMUM  DIVERGENCE  RATES  FOR  THREE-BLADED  ROTOR 


Helix  (rad) 


1.75 
2.00 

2.25 

2.50 
3.00 

3.50 

3.75 
4.00 

4.50 
5.00 

5.25 

5.50 
6.00 
0.50 


Max  Divergence  Rate 


113.97 

100.98 
78.83 
94.00 

115.57 
123.71 

116.57 
96.52 
79.65 

96.85 
115.37 

99.41 

79.28 

90.88 

114.81 

95.85 
79.47 
89.97 

108.00 

45.38 

41.62 

54.67 

43.79 

50.17 

38.36 

97.19 

93.55 

120.01 

95.78 

114.16 

89.40 
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TABLE  VIII.  MAXIMUM  DIVERGENCE  RATES  FOR  FOUR-BLADED  ROTOR 


Pitch  of 
Helix  (rad) 


Max  Divergence  Rate 


209.37 

194.22 

168.98 

135.16 

165.99 
191.69 
207.28 
212.87 
193.32 

138.17 
^  186.24 

208.16 

190.05 

139.94 

179.51 

201.23 
141.21 

193.23 
87.25 
61.20 
84.72 
94.52 
63.58 
90.31 
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TABLE  IX.  MAXIMUM  DIVERGENCE  RATES 

FOR  FIVE-BLADED  ROTOR 

Pitch  of 

Helix  (red) 

Wave 

Number 

Core 

Size 

Max  Divergence  Rate 

0.10 

0.25 

0.10 

332.56 

0.10 

o.so 

0.10 

315.78 

0.75 

0.10 

288.09 

0.10 

1.00 

0.10 

250.72 

0.10 

1.25 

0.10 

219.68 

0.10 

0.10 

1.50 

1.75 

jjjgj* 

252.93 

282.58 

0.10 

2.00 

0.10 

307.43 

0.10 

2.50 

0.10 

340.74 

0.10 

3.00 

0,10 

311.15 

0.10 

3.50 

0.10 

256.42 

0.10 

3.75 

0.10 

220.27 

0.10 

4.00 

0.10 

242.70 

0.10 

4.50 

0.10 

302.72 

0.10 

5.00 

0.10 

327.21 

0.10 

5.50 

0.10 

307.83 

0.10 

6.00 

0.10 

253.02 

0.10 

6.25 

0.10 

220.94 

0.10 

7.00 

0.10 

290.84 

0.10 

7.50 

0.10 

316.34 

0.10 

8.00 

0.10 

298.76 

TABLE  X.  MAXIMUM  DIVERGENCE  RATES  FOR  SIX-BLADED  ROTOR 


Pitch  of 
Helix  (rad) 

Wave 

Number 

Core 

Size 

Max  Divergence  Rate 

0.10 

0.25 

■Hflj 

483.43 

0.10 

0.50 

465.33 

0.10 

0.75 

435.62 

0.10 

1.00 

HQ 

395.08 

0.10 

1.25 

358.51 

0.10 

1.50 

0.10 

318.15 

0.10 

1.75 

0.10 

353.64 

0.10 

2.00 

0.10 

390.09 

0.10 

2.50 

0.10 

459.22 

0.10 

3.00 

0.10 

485.01 

0.10 

3.50 

0.10 

462.12 

0.10 

4.00 

0.10 

396.97 

0.10 

4.50 

0.10 

318.72 

0.10 

5.00 

0.10 

381.72 

0.10 

5.50 

0.10 

445.94 

0.10 

6.00 

0.10 

472.46 

0.10 

7.00 

0.10 

394.10 

0.10 

8.00 

0.10 

370.54 
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In  short,  peak  values  of  divergence  rate  occur  at  wave  numbers, 
co,  equal  to  integer  multiples  times  one-half  the  number  of 
blades,  and  in  all  cases  the  magnitude  of  the  peak  divergence 
rate  decreases  with  wave  number. 

Enveloping  curves  can  be  drawn  through  maximum  and  minimum  di¬ 
vergence  rate  points  in  all  cases.  There  is  some  absolute  max¬ 
imum  divergence  rate  which  occurs  at  the  maximum  of  the  envel¬ 
oping  curve  through  the  maximum  divergence  rate  points.  For 
two-,  four-  and  six-bladed  rotors  this  absolute  maximum  rate 
occurs  at  wave  number  0.0,  i.e.,  for  the  dilatational  mode. 

For  three-  and  five-bladed  rotors  the  dilatational  mode  is 
still  highly  unstable,  but  the  absolute  maximum  divergence  rate 
occurs  at  wave  numbers  l.S  and  2.5  respectively,  i.e.,  where 
the  wave  number  is  one-half  the  number  of  blades  (  to  »  A/2). 

To  assess  the  effect  of  helix  pitch  angle,  some  variations  in 
this  parameter  were  run  for  two-,  three-, and  four-bladed  ro¬ 
tors.  For  the  higher  pitch  angles  the  divergence  rate  drops 
considerably,  even  though  the  characteristic  variations  with 
wave  number  remain.  For  example,  an  increase  in  pitch  of  helix 
from  0.10  to  0.1S  reduces  the  divergence  rate  by  a  factor  of 
almost  2.  Note  that  the  curves  for  pitch  *  0.15  must  be  read 
on  the  right  scales,  as  also  pointed  out  in  Figures  15  through  17. 
Physically,  this  behavior  can  be  explained  in  terms  of  inter¬ 
action  between  the  neighboring  coils  of  the  helices.  For  high¬ 
er  pitch  angles  the  neighboring  helix  coils  are  farther  apart, 
and  distortions  that  tena  to  cause  neighboring  coils  to  ap¬ 
proach  each  other  therefore  result  in  lesser  interactions.  Ab¬ 
solute  maximum  divergence  rates  for  different  rotors  are  plotted 
against  pitch  in  Figure  22  and  are  tabulated  in  Table  XI. 

The  nondimensional  divergence  rate  of  multiple-bladed  rotors 
would  seem  to  increase  with  the  number  of  blades  for  the  same 
sort  of  reason  that  it  increases  with  decreasing  pitch;  namely, 
there  is  a  decreased  separation  between  successive  lines  of 
vorticity.  Figure  22  shows  that  for  a  helix  pitch  angl£  of 
.05,  for  example,  a  two-bladed  rotor  has  a  maximum  divergence 
rate  of  200  and  a  four-bladed  rotor  has  a  maximum  divergence 
rate  of  850.  It  must  be  noted  that  these  rates  are  nondimen- 
sionalized  by  total  vortex  strength  r  ,  as  shown  in  the  equa¬ 
tion  for  A  on  pages  27  and  37.  If  the  sum  of  the  blades  in 
each  rotor  produced  the  same  lift  (which  is  the  case  in  compar¬ 
ing  alternate  designs  for  the  same  aircraft),  the  vortex 
strength,  V  ,  in  the  two-bladed  case  would  be  twice  that,  of  the 
four-bladed  case,  so  that  the  dimensional  rates  would  then 
compare  as  400  to  850.  Thus,  in  terms  of  the  stability  of  the 
vortices  in  the  far  wake  for  rotors  of  equal  lifting  capacity, 
the  fewer  blades  the  better. 

The  effect  of  core  size  on  stability  is  not  so  noticeable  for 
multiple  rotors.  When  the  core  size  was  increased  from  0.1  to 
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TABLE  XI 


ABSOLUTE  MAXIMUM  DIVERGENCE  RATES  FOR 
MULTIBLADED  ROTORS  OF  VARYING  PITCH 


Pitch  of  Number  of 
Helix  (rad)  Blades 


Wave 

Number 


Absolute  Maximum 
Divergence  Rates 


0.  05 
0.10 
0.15 
0.20 


0.0 


199.89 

50.05 

22.34 

12.65 


0.05 

0.10 

0.15 

0.20 


1.5 


492.11 

123.71 

54.67 

30.47 


0.05 

0.10 

0.15 

0.20 


0.0 


853.90 

214.57 

96.17 

54.70 


0.05 
0.10 
0.15 
0. 20 


1363.00 

340.74 

150.64 

84.26 


0.05 

0.10 

0.15 

0.20 


1940.50 

489.65 

219.78 

125.26 


0.33,  the  amplification  rate  dropped  only  about  21  for  the  wave 
numbers  with  highest  divergence  rates  in  all  cases  where  r\  ■  1. 
Core  size  does  appear  to  be  an  important  influence  on  stability 
of  a  single  helix,  i.e.,  where  n.  ■  1,  as  shown  in  Figure  14.. 
Its  reduced  effect  on  the  stability  of  multiple-bladed  rotors 
does  not  seem  unreasonable,  since  the  mutual  inductance  between 
different  helices  is  a  relatively  important  influence  on  the 
stability  of  interdigitated  vortex  systems,  whereas  variations 
in  core  size  are  felt  most  in  self-induction  effects. 


4.8 


COMMENTS  ON  ABSOLUTE  MAXIMUM  WAKE  DIVERGENCE  RATES  FOR 

MULTIBLAPEP  ROTORS  VITH~PTO  XKP  my  NUMBERS  6 F  TTaEEs 


This  section  deals  with  a  qualitative  physical  explanation  of 
absolute  maximum  divergence  rates  for  multibladed  rotors. 

First,  as  seen  in  Figures  15  through  19,  these  absolute  maxima 
are  associated  with  the  0  and  ■  n/2  modes  for  rotors 

with  even  and  odd  numbers  of  blades  respectively.  Comments  in 
most  cases  will  be  restricted  to  axial  perturbation  displace¬ 
ments  since  the  effect  on  induced  velocities  of  such  motions 
is  one  order  of  magnitude  higher  than  either  radial  or  circum¬ 
ferential  perturbation  displacements  insofar  as  the  mutual 
inductance  between  neighboring  coils  is  concerned.  This  can 
be  explained  more  explicitly  as  follows: 

The  Biot-Savart  law  reveals  that  the  maximum  interaction  occurs 
between  those  vortex  elements  which  are  nearest  to  each  other. 
For  helix  pitch  angles  less  than  0.20,  the  parallel  elements 
which  are  closest  to  each  other  are  on  two  neighboring  coils 
at  the  same  circumferential  position,  i.e.,  located  on  the  same 
line  parallel  to  the  axial  direction.  This  axial  distance 
shall  be  denoted  byd.  A  radial  perturbation  of  in  the  radial 
direction  will  replace  d.  by  ,  and  the  integrand  in  the 

Biot-Savart  law  will  be  modified  from  d**  to  dT*  C‘<- 3</d  +  -  j, 

contributing  a  second -order  term  to  the  perturbation  integrand. 
However, an  axial  perturbation  of  £  will  modify  the  minimum 
interaction  distance  d.  to  td.- t>,  contributing  Cd-t)'4  in  the 
perturbation  integrand.  This  corresponds  to  a  change  from 

d'3  to  ct  [|+a4-^+ .  ]  and  is  a  first-order  contribution. 

Therefore,  the  axial  perturbation  is  the  dominant  factor  in  the 
interactions  between  the  nearest  parallel  vortex  elements.  It 
follows  that  the  divergence  rate  of  a  given  deformation  mode 
will  be  greater  than  another,  if  the  normalized  axial  displace¬ 
ments  which  reduce  the  separation  between  vortex  elements  adja¬ 
cent  to  their  counterparts  on  neighboring  coils  take  place 
over  more  of  the  arc  lengths  of  the  helix  than  the  other  defor¬ 
mation  mode. 
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For  the  purpose  of  evaluating  this  intuitive  notion  pertaining 
to  the  absolute  maximum  wake  instability  of  multibladed  hover¬ 
ing  rotors,  an  average  distortion  factor  will  be  defined. 

This  ’’average  relative  distortion”  will  be  calculated  for  two- 
and  three-bladed  rotors  for  axial  perturbations  at  the  wave 
numbers  where  maximum  instabilities  occur.  Subsequently,  it 
will  be  used  to  explain  why  absolute  maximum  instability 
occurs  at  wave  numbers  0.0  for  rotors  with  an  even  number  of 
blades,  and  at  wave  numbers  1.5  and  2.5  for  rotors  with  an  odd 
number  of  blades,  i.e.,  three-  and  five-bladed  rotors, 
respectively. 

For  the  case  of  the  three-bladed  rotor,  maximum  instabilities 
occur  at  wave  numbers  0.0  and  1.5.  Eigenfunctions  for  the 
most  unstable  eigenvalues  at  these  wave  numbers  are  tabulated 
in  Table  XII, and  the  elevation  views  of  the  perturbed  vortex 
system  for  only  axial  perturbations  are  sketched  in  Figures 
23  and  24. 


For  a  wave  number  of  0.0,  eigenfunctions  for  helices  1,  2  and 
3  are  a  a  a 

it,  «  +  l  ;  Sii  =  O  Sij*  -  l 

The  "average  relative  distortion"  is  defined  as  a  line  integral 
of  the  reduction  in  the  axial  displacement  between  vortex 
elements  at  the  same  circumferential  position  on  neighboring 
coils  evaluated  around  the  azimuth. 


A  M  J 

Thus,  the  average  relative  distortion  between  the  1  and  3 

‘  ‘  '  '  ‘  .71 


helices  for  period  ir is 


(-*  f 

Jo  j 


ale  r  zn 

a 


The  average  relative  distortion  between  the  3rt*  and  2n<*  helices 
for  period  rr  is 


1 


rr 

oie  ■=  tt 


and  finally,  the  average  relative  distortion  between  the  2 


nd 


and  lsthelices  for  period  it  is 


dO 


The  total  average  relative  distortion  factor,  therefore,  is 
4tt  -  12.57. 

For  a  wave  number  of  1.5,  eigenfunctions  for  helices  1,  2  and 

t  arp  a  a  A 

S6|e  -*-1  J  St*  =  -  0.25  j  =  V  1 
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TABLE  XII.  EIGENFUNCTIONS  FOR  THE  MOST  UNSTABLE  EIGENVALUES 


OF  A  THREE-BLADED  ECTOR  AT  WAVE  NUMBERS  0.0 
AND  1.5,  PITCH  -  0.1  AND  CORE  SIZE  -  0.1 


Eigenfunction  for 

Wave  Number  0.0 

Eigenfunction  for 

Wave  Number  1 . 5 

Label 

Real 

Imaginary 

_ 

Label 

Real 

Imaginary 

A 

6r, 

-1.000 

BH 

|  B 

-0.076 

B 

St, 

♦0.996 

■Bp  ■ 

A 

S*. 

+0.996 

+0.004 

A- 

1*2. 

0.000 

+0.210 

0.000 

0.000 

0.000 

rtA; 

+0.016 

+0.002 

0.000 

0.000 

-0.248 

+0.000 

A 

1.000 

C.000 

-0.981 

-0.005 

0.076 

0.000 

t4s 

-0.073 

-0.012 

A 

-0.996 

0.000 

+0.996 

+0.004 
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Referring  to  Figure  24,  the  following  distortion  factors  can 
be  calculated. 

The  average  relative  distortion  between  1st  and  3rd  helices 
for  period  tt  is 

C«'t  rW/* 

1 1 ««'«i (!•<#>  I  f  coiose)l]<*e  4-iJ  ClCoso  «a>l  4-  J^O 

IT/3  f^3 

+  Cl  cot  0-70)1  -  is.v,(«-se>n<io^Jcitosc..«e)l-i6ivuis#)|]<t© 

TT/ j  f» 71/3 

+  Jcicosoee)l  +|  s.y\ct-60)l]oie  vi  Jc.lcos*).se)l-HSin(t-5©>l  jeie 

■  4-  +  4.C 

Similarly,  the  average  relative  distortion  between  3rd  and  2nd 
helices  for  period  ir  is 

r  r  ^S  T  r- 

J CK&s  (1.5031  +  0.25 1  StV)  o-se  >1 3°^©  * JC  |Co*(I'5©)1  +  0.2S|Si,«(i-MHJ4® 

r  %  r 

♦  J  DS''n('Se)1  -0.2S  l  cot  (i.5e)|  ]el6v  l  J  ft6inCi.5e)|-0.2.S|C»li  S©)j^le 

rn/*  fTT-/» 

+  J  tlCoS(i-5e)|4.0.25lS.V)(i.5e)Ocle  +■  lJCICo*0.5e>l+0.25lSiVKI*S8j3eU 

*  3  2S  +  5-lS  i 

Finally,  the  average  relative  distortion  between  2nd  and 
1st  helices  for  period  TT  is 

Jco.2S|6tVi((.5S)l+  I  5iV>Cl.5e)|]cL©  +■  i  JCo.aS  I Sir>ci  69)1+1 

r  ^ 

icoitisejf  +  o.it|Cos(i.se)l]i0  +  i J C.I 

J0  0 

rr/i  r 

+  J[o.2ff  IS.nOS«)»  +  I  S,V»OS«l>t®  +  iJ^0.15|3iv»Cl.5©)|  +!&!«(!  50)1]^© 

«  3-?S,4.3.>Si 
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Thus .  ths  total  average  ralativa  distortion  factor  is  11  ♦  Hi, 
and  its  magnitude  is  approximately  15.5. 

Consider  the  case  of  a  two-bladed  rotor;  maximum  instabilities 
occur  at  wave  numbers  0.0  and  1.0.  Eigenfunctions  for  the 
most  unstable  eigenvalues  at  these  wave  numbers  are  tabulated 
in  Table  XIII. 

For  a  wave  number  of  0.0,  the  axial  distortion  amplitudes  for 
helices  1  and  2  are 

St,  *  -  i.oo  j  loo 

The  average  relative  distortion  factor  over  a  period  it  is 
Cn  f * 

Cle  4-  de  *  2TT  Sf  6.26 

For  a  wave  number  of  1.0,  the  axial  distortion  amplitudes  for 
helices  1  and  2  are 

6i, -o  tu  )  0  612 

Referring  to  Figure  21,  the  average  relative  distortion  factor 
over  a  period  IT  is  _ 

r  r’W 

2(011 2)  icotcexde  4-  2(osi2)  I  coS(»;l<4e 
Jo  Jo 

+  i[2.coeia)Ji*T'n<e;i«ie  4-  2(oau)J)si»i(e)» 

■  3.1R-*  V«l* 

Thus  the  magnitude  of  the  average  relative  distortion  factor 
is  approximately  4.60. 

Results  of  the  above  calculations  are  tabulated  in  Table  XIV. 
This  shows  that  the  relative  axial  displacements  do,  in  fact, 
determine  the  stability  and  that  those  modes  which  have 
the  highest  average  reduction  of  axial  separation  between 
adjacent  coils  around  the  circumference  of  the  helix  will  be  the 
most  highly  divergent. 
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TABLE  XIII.  EIGENFUNCTIONS  FOR  THE  MOST  UNSTABLE  EIGENVALUES 
OF  A  TWO-BLADED  ROTOR  AT  WAVE  NUMBERS  0.0 
AND  1.0,  PITCH  *0.1  AND  CORE  SIZE  -  0.1 


Eigenfunction  for 
Wave  Number  0.0 


Eigenfunction  for 
Wave  Number  1.0 


Label 


Imaginary 


Label 


Imaginary 


+0.998 

+0.065 

-0.998 

-0.998 

-0.065 

+0.998 


-0.002 

-0.002 

-0.000 

+0.002 


+0.783 

-0.784 

-0.043 

+0.783 


-0.215 

-0.034 

+0.212 
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TABLE  XIV.  COMPARISON  OF  MAXIMUM  INSTABILITY  FOR 
TWO-  AND  THREE -BLADED  ROTORS 


Number 

of 

Blades 

Wave 

Numbers 

Amplitude  of 
Average  Dis¬ 
tortion  Factor 

Max  Di¬ 
vergence 
Rate 

Conclusion 

2 

0.0 

6.28 

50.05 

Perturbations 
for  the  wave 
number  0.0  are 
more  unstable 

1.0 

4.60 

49.80 

3 

0.0 

12.57 

118.42 

Perturbations 
for  the  wave 
number  1.5  art- 
more  unstable 

1.5 

15.50 

123.71 
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CHAPTER  S 


CONCLUSIONS 


1.  The  small-perturbation  stability  analysis  of  a  doubly  infi¬ 
nite  array  of  interdigitated  right  circular  helical  vortices 
formulated  here  shows  reasonable  agreement  with  the  recent 
work  of  Widnall  [Reference  16]  where  vortex  core  to  helix 
diameter  ratios  are  0.1  or  less  and  when  the  number  of  blades 
is  set  equal  to  unity.  The  larger  differences  which  appear  in 
a  comparison  with  the  Reference  16  results  at  larger  core 
ratios  (e.g.,  0.33)  are  attributed  to  differences  in  the 
methods  for  eliminating  the  self-induction  singularity.  While 
the  method  works  best  for  small  pitch  angles,  it  has  been 
shown  to  be  the  equivalent  of  the  formulations  developed  by 
S.C.  Crow  and  Parks  [References  17  and  20]  dealing  with  the 
infinite  straight  vortices  such  as  trail  from  the  tips  of 
fixed-wing  aircraft,  when  the  pitch  angle  is  made  to  approach 
infinity  for  a  two-bladed  rotor. 

2.  A  continuum  of  instability  modes  has  been  found  associated 
with  all  values  of  wave  numbers;  only  modes  with  wave  numbers 
0  and  1  are  so  much  as  neutrally  stable,  and  only  for  the  case 
of  a  single  helix.  The  most  unstable  modes  involve  the  most 
axial  motion  of  adjacent  vortex  segments  relative  to  each  other. 
By  "adjacent  segments"  is  meant  vortex  segments  above  and  below 
each  other  (i.e.,  at  the  same  azimuthal  location)  on  adjacent 
coils  of  the  same  or  neighboring  helices.  Furthermore,  the 
larger  the  percentage  of  the  helical  arc  length  involved  in 
such  motion,  the  more  rapidly  the  distortion  will  diverge. 
Maximum  divergence  rates  in  the  unstable  modes  increase  as  the 
helix  pitch  decreases,  increase  as  the  number  of  helices 
increases,  and  decrease  as  the  number  of  cycles  of  deformations 
in  one  turn  of  the  helix  (i.e.,  wave  number)  increases. 

3.  The  effect  of  increasing  helix  filament  core  diameter  is  to 
make  the  analysis  more  sensitive  to  the  means  by  which  the  sing¬ 
ularity  is  eliminated.  Increasing  the  number  of  blades  appears 
to  make  the  core  diameter  less  influencial  in  determining  the 
rates  of  divergence.  Although  it  has  not  been  demonstrated 
here,  one  would  expect  that  reduced  pitch  would  also  reduce  the 
dependence  of  the  calculation  on  the  vortex  core  size.  Finally, 
a  larger  core  diameter  ratio  also  appears  to  reduce  the 
difference  between  peak  and  minimum  divergence  rates  as  a 
function  of  wave  number. 
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CHAPTER  6 


RECOMMENDATIONS 


In  view  of  the  instability  of  helical  vortex  arrays  observed 
(a)  in  tests  and  (b)  in  iterative  distorted  wake  geometry 
calculations  and  which  now  have  been  confirmed  in  this  study, 
it  is  recommended  that  the  stability  of  a  semi-infinite  vortex 
array,  together  with  the  bound  vortices  which  generate  it,  be 
studied  in  a  manner  similar  to  the  analyses  reported  here. 

The  purpose  of  such  a  study  would  be  twofold:  first,  to 
examine  the  effect  of  the  bound  vortex  on  helical  vortex 
stability  and  the  associted  unstable  divergence  rates,  and  . 
second,  to  develop  analytical  and  realistic  perturbation 
quantities  to  be  used  as  forcing  functions  for  the  unstable 
modes  of  the  semi- infinite  vortex  array  typical  of  lifting 
rotors.  In  such  a  "forced"  calculation,  one  would  expect 
to  predict  the  actual  distorted  geometry  of  the  wake  beneath 
lifting  rotors  which  must  be  known  to  perform  more  accurate 
airloads  analyses,  for  all  purposes. 
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Figure  1.  Undistorted  p^and*?1  Vortex  Helices  of  Radius  r.  Pitch  ft 
and  Vortex  Core  Diameter 
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Figure  2.  Plot  of  Divergence  Rate  Against  Phase  Difference 
for  a  Two-Bladed  Rotor. 


Figure  3.  Major  Steps  in  the  Computer  Program 

72 


74 


025 


Against  * ; 

-kx.l<u*»)x*«.e>} **)  Against 


Figure  7.  Radial  Perturbations  in  Wave  Number  ■  0  Mode 
(Dilatational) . 
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Figure  8.  Radial  Perturbations  in  Wave  Number  ■  l  Mode. 
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Figure  9.  Axial  Perturbations  in  Wave  Number  *  1  Mode. 
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Figure  11.  Axial  Perturbations  in  Nave  Number  "  2  Mode. 
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Figure  12.  Radial  Perturbations  in  Wave  Nuaber  ■  1/2  Mode, 
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Figure  14.  Plot  of  Maximum  Divergence  Rate  vs  Wave 
Number  for  Single  Helix. 
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Figure  15.  Flot  of  Maximum  Divergence  Rate  vs  Wave 
Number  for  Two*Bladed  Rotor. 
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Figure  16.  Plot  of  Maximum  Divergence  Rate  vs  Wave 
Number  for  Three- Bladed  Rotor. 
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Figure  IT.  plot  of  Maximum  Divergence  Rate  vs  Wave 
Number  for  Four-Bladed  Rotor. 
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Figure  18.  Flot  of  Maximum  Divergence  Rate  vs 
Number  for  Five-Bladed  Rotor. 
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Figure  19.  Plot  of  Maximus  Divergence  Rate  vs 
Number  for  Six~Bladed  Rotor. 


80 

Wave 


88 


Two-Bladed  Rotor 


VO 

mT 


V2* 


Two-Bladed  Rotor 


0.0 


Figure 


005  OIO  015  040 

Pitch  of  htlix  (Rodion*)—— 


22.  Absolute  Maximum  Divergence  Rates  for 
Multibladed  Rotors  of  Varying  Pitch. 
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Figure  23.  Real  Component  of  Axial  Perturbations  for  a 
Three-Bladed  Rotor  at  Wave  Number  0.0. 


92 


LITERA  .'URL  CITED 


1.  Glauert,  II.,  AIRPLANE  PROPFLLF.RS ,  Aerodynamic  Theory, 
edited  by  K.  F.  Durand,  Dover  Publications,  New  York,  1935, 
Chapter  X,  Division  L,  Vol.  IV. 

2.  Shapiro,  Jacob,  PRINCIPLES  OF  HFLICOTTER  ENGINEERING,  New 
York,  McGraw-Hill  Book  Co.,  Inc.,  1956,  pp.  70-73. 

3.  Gessow,  A. .and  Myers,  G.  C.  Jr.,  AERODYNAMICS  OF  THE  HELI¬ 
COPTER,  New  York,  The  MacMillan  Company,  1952,  pp.  46-47. 

4.  Mangier,  K.  K. ,  CALCULATION  OF  INDUCED  VELOCITY  FIELD  OF  A 
ROTOR,  Royal  Aircraft  Establ i shmcnt  Report,  No.  Aero.  2247, 
1948. 

5.  Loewy,  R.,and  Joglekar,  M. ,  AN  ACTUATOR  DISC  ANALYSIS  OF 
HELICOPTER  WAKE  GEOMETRY  AND  THE  CORRESPONDING  BLADE  RESPONSE, 
University  of  Rochester;  USAAVLABS  lechnical  Report  69-66, 

U.  S.  Army  Aviation  Materiel  Laboratories,  Fort  Eustis, 
Virginia,  December  1970,  AD  891981. 

6.  Goldstein,  S.,  ON  THE  VORTEX  THEORY  OF  SCREW  PROPELLERS, 
Proceedings  -  Royal  Society  of  London,  Vol.  A123,  No.  A792, 

7.  Jenny,  D. ,  Olson,  J.,  and  Landgrebe,  A.,  A  REASSESSMENT  OF 
ROTOR  HOVERING  PERFORMANCE  PREDICTION  METHOES,  Journal  of  the 
American  Helicopter  Society,  Vol.  13,  No.  2,  April  1968, 

pp.  1-Z6. 

8.  Miller,  R.  H.,  ON  THE  COMPUTATION  OF  AIRLOADS  ACTING  ON 
ROTOR  BLADES  IN  FORWARD  FLIGHT,  Journal  of  the  American  Heli¬ 
copter  Society,  Vol.  7,  No.  2,  April  196Z. 

9.  DuWaldt ,  F.  A.,  Pitiali,  R.  A.,  A  METHOD  FOR  COMPUTING 
ROTARY  WING  AIRLOAD  DISTRIBUTIONS  IN  FORWARD  FLIGHT,  U.  S. 

Army  TRECOM  Report  No.  TCREC  TR  62-44,  November  1962. 

10.  Tarrarine,  S.,  EXPERIMENTAL  AND  THEORETICAL  STUDY  OF  LOCAL 
INDUCED  VELOCITIES  OVER  A  ROTOR  DISC  FOR  ANALYTICAL  EVALUATION 
OF  THE  PRIMARY  LOADS  ACTING  ON  HELICOPTER  ROTOR  BLADES,  Euro¬ 
pean  Research  Office,  U.  S.  Army  Contract,  October  1960. 

11.  Landgrebe,  A.  J.,  AN  ANALYTICAL  AND  EXPERIMENTAL  INVESTI¬ 
GATION  OF  HELICOPTER  ROTOR  HOVER  PERFORMANCE  AND  WAKF.  GEOMETRY 
CHARACTERISTICS,  United  Aircraft  Corporation;  USAAMRDL  Tech¬ 
nical  Report  71-24,  U.  S.  Air  Mobility  Research  and  Develop¬ 
ment  Laboratory,  Fort  Eustis,  Virginia,  June  1971,  AD  728835. 


94 


12.  Crimi,  P.,  PREDICTION  OF  ROTOR  WAKE  FLOKS,  presented  at 
CAL/USAAVLABS  symposium,  Buffalo,  New  York,  June  1966. 

13.  Sadler,  Gene  S.,  MAIN'  ROTOR  FRFF  KAKF  GEOMETRY  EFFECTS  ON 
BLADE  AIR  LOADS  AND  RFSPONSr  TOR  HELICOPTERS  IN  STEADY 
MANEUVERS,  Rochester  Applied  Science  Associates,  Inc.,  NASA 
CR-2110,  National  Aeronautics  and  Space  Administration, 
Washington,  D.  C. ,  September  1972. 

14.  Kfrmin,  Theodore  von,  AERODYNAMICS,  Ithaca,  N.  Y.,  Cornell 
University  Press,  19S4,  Chapter  II,  pp.  31-59. 

15.  Levy,  M.  A.,  and  Forsdyke,  A.  G.,  THE  STEADY  MOTION  AND 
STABILITY  OF  A  HELICAL  VORTEX.  Royal  Society  -  Proceedings. 

Vol.  120,  No.  A786 ,  October  1,  1928,  pp.  670-690. 

16.  Kidnall,  S.  E.,  THE  STABILITY  OF  A  HELICAL  VORTEX  FILAMENT, 
Journal  of  Fluid  Mechanics.  Vol.  54,  Part  4,  1972,  pp.  641-663. 

17.  Crow,  S.  C.,  STABILITY  THEORY  FOR  A  PAIR  OF  TRAILING 
VORTICES,  American  Institute  of  Aeronautics  and  Astronautics 
Journal ,  Vol.  S,  Mo.  12,  Dec.  197CT, pp.  "2l72-2l7l9. 

18.  Hama,  F.  R. ,  PROGRESSIVE  DEFORMATION  OF  A  CURVED  VORTEX 
FILAMENT  BY  ITS  OWN  INDUCTION,  The  Physics  of  Fluids.  Vol.  S, 
No.  10,  October  1962,  pp.  11S6-116?.- 

19.  Hama,  F.  R. ,  PROGRESSIVE  DEFORMATION  OF  A  PERTURBED  LINE 
VORTEX  FILAMENT.  The  Physics  of  Fluids.  Vol.  6,  No.  4,  April 
1963,  pp.  526-534. 

20.  Parks,  P.  C.,  A  NEW  LOOK  AT  THE  DYNAMICS  OF  VORTICES  WITH 
TINITE  CORES:  Proceedings  of  Symposium  on  Aircraft  Wake  Tur¬ 
bulence  held  in  Seattle,  Washington,  September  1970  (Included 
in  the  book  AIRCRAFT  WAKE  TURBULENCE  AND  ITS  DETECTION, 
editors,  Olson,  J.  H.,  Goldberg,  A., and  Rogers,  M. ,  New  York, 
Plenum  Press,  1971,  pp.  355-388). 

21.  Batchelor,  G.  K.,  AN  INTRODUCTION  TO  FLUID  DYNAMICS, 
Cambridge  University  Press,  Cambridge,  England,  1967,  p.  264. 

22.  Hildebrand,  F.  B.,  INTRODUCTION  TO  NUMERICAL  ANALYSIS, 

New  York,  McGraw-Hill  Book  Company,  1956,  pp.  73,  141,  146. 

23.  Krylov,  V.  I.,  APPROXIMATE  CALCULATION  OT  INTEGRALS,  Ne* 
York,  The  Macmillan  Company,  1962,  pp.  338-339. 

24.  Ralston,  A.,  A  FIRST  COURSE  IN  NUMERICAL  ANALYSIS,  New 
York,  McGraw-Hill  Book  Company,  1965,  pp.  516-519. 


95 


25.  Acton,  F.  S.,  NUMERICAL  METHODS  THAT  WORK,  New  York, 
Harper  and  Row,  1970,  p.  347. 

26.  Payne,  P.  R.,  HELICOPTER  DYNAMICS  AND  AERODYNAMICS,  New 
York,  The  Macmillan  Company,  1959,  p.  32. 

27.  Thomson,  W.,  ON  THE  VIBRATIONS  OF  A  COLUMNAR  VORTEX, 

Philosophical  Magazine.  Serial  5,  Volume  10,  No.  61,  Sept. 

1580,  pp.155-168. - 

28.  Widnall,  S.  E.,  and  Sullivan,  J.,  ON  THE  STABILITY  OF 

VORTEX  RINGS,  Proceedings  -  Royal  Society  of  London.  Vol. 
A332,  1973,  pp.  335-3537 - - - 

29.  Betchov,  R.,  ON  THE  CURVATURE  AND  TORSION  OF  AN  ISOLATED 
VORTEX  FILAMENT.  Journal  of  Fluid  Mechanics.  Vol.  22.  Part  3, 
1965,  pp.  471-479” 

30.  Spreiter,  J.  R.,and  Sacks,  A.  H.,  THE  ROLLING  UP  OF  THE 
TRAILING  VORTEX  SHEET  AND  ITS  EFFECT  ON  THE  DOWNWASH  BEHIND 
WINGS,  Journal  of  the  Aeronautical  Sciences,  Vol.  18,  No.  1, 
Jan.  1 


96 


APPENDIX  I 

LIMITS  FOR  S.  C.  CROW'S  CASE 


The  equivalence  of  the  special  case  of  the  present  work  with 
that  of  Crow  and  Parks  can  be  seen  by  noting  that  the  direc¬ 
tion  z  in  the  present  study  corresponds  to  the  axial  direction 
of  the  aircraft  trailing  vortex,  and  that  the  directions  r  and 
*•  of  the  present  study  in  the  limit  will  correspond  to  the  y 
and  a  directions,  respectively,  in  S.  C.  Crow's  theory.  All  the 
nondinensional  parameters  in  Crow's  work  must  be  defined  in 
terms  of  the  parameters  of  the  present  study  to  demonstrate  the 
correspondence.  The  r  perturbation  equation  in  the  present 
study  with  above-mentioned  restrictions  can  be  written  as 

*  J^1,  [  1  s;*  i 

Mai 


-  Cos  («•»  p  4  Cos  (OXwf  ♦  S.W  (  vVm^JCoS 

4-  to  K«.r  )  Co%‘QX-*~}  -  Co* 

♦  S!«  (*-*♦  qw^^a.-ncox^  \  ♦  fcv'adu  <•  )Co*i«o  *<-*> 


4-  X*«*  Si<*  (  +  iy^>Cot(4)X^ 

4-  ICX^  Cose  1  * 

54.  C*»  (.*->♦  6M CO Km^  4-  X***  &V»CX«,^4.»P*.^  S,«  <0*«vtp 

*  A 

-  «•>  Mwf  C«k(,x<m^  4  4W*  JCcS03K*»>y  )  4  *  Vv  C®*  t  Hw*y  4  MR-*) 


-  Xm»»  c  *  vV-*  )  J  3  eLx—p 


VV  -  XlOf 


CoS  (.It-* 4.  •)}] 


where  Jmp  and  are  defined  in  Section  2.6. 

First  the  variables  of  the  present  theory  will  be  transformed 
into  corresponding  variables  used  in  S.  C.  Crow's  analysis, 
and  then  the  limit  will  be  taken  as  helix  pitch  tends  to 
infinity.  All  of  S.  C.  Crow's  parameters  will  carry  the 
subscript  's'  to  avoid  confusion. 

From  the  definition  of  axial  coordinate,  Rrx»Xe. 

« t  4  «■  (*9  X— 4 

The  perturbations  in  the  present  study  are  of  the  type  e 
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and  if  the  wave  number  in  S.  C.  Crow's  theory  is  expressed  as 
fcj  »  oz.  ,  the  form  of  the  perturbations  is  transformed  to 


nr 


and  the  nondimens ional  wave  number  becomes 


«  lu  tr  :  ktb  . 

Applying  the  transformation  described  above,  the  self- indue* 
tion  term  can  be  written  as 

«*r,  ,  £  [il-  S;Wj£) 

o 

*R  £V*  hr'ti  6»  6'v'(.ft)C"t«r‘) 

fc*1  kr‘‘5.  c-  ^(sJ)  -  g*  so  (.&)]  ^ 

-  «o 

-  5/2 

....  plus  terms  involving  Jjj  •' 

Now,  *>  m  x.  ce-»o  as  helix  pitch  -»eo  ,  and  one  can  approxi¬ 
mate 


Cot 


**  (»)  -  s? 


Neglecting  terms  and  noting  that  Ri,  it  ^ 

follows  that  Cat'-  ar'Co.(&>  5^*V  «!]  •  j  ‘ 

•nd  w  tr,  .  ]■*». 

Similarly,  mutual -inductance  terms  can  be  written  as 
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*  ii  It,***  [Co*c5>^n)  CosifciX*)  ♦  ♦  rr)  Cos(fe*x».) 

Jo 

•*■  *s*s  Cot(-S^'  vTT)  S,V*  ***  *0  Tt^ 

4  s  c-costigt  +*)  -  *»  s>-  c^r^)3^‘ 
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Taking  the  limits  again  yields 
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Writing  s£.5,t 


»**$,***.  ¥*»  and 


and  remem¬ 


bering  that  fcr  is  positive  outward,  while  y.  and  y,.  are  paral¬ 
lel  (thereby  changing  sign  of  either  the  self-  or  mutual- 
inductance  terns), 
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This  is  the  same  equation  that  results  when  Parks'  modification 
is  applied  to  Crow's  theory.  The  one  remaining  equation 

relating  £,*  with  y(%  and  y»*  can  be  obtained  by  following  the 
sane  process  for  y&$(  perturbation  equation. 
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APPENDIX  II 

FORTRAN  IV  SOURCE  LISTING 


Main  Program . 101 

Subroutine  ALLMAT  .  112 

Subroutine  QG9 . 121 

Subroutine  INTG . 123 

Function  Subprogram  DSIMP . 130 
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